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INTRODUCTION 

River-basin development has become one of the largest classes of 

public enterprise in the United States. Multiple-purpose river-basin 
programs may involve power, irrigation, flood control, pollution 
control, navigation, municipal and industrial water supply, rccrca- 
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tion, fish and wildlife, and the conservation of soil and water on 
watershed lands. 

Almost every kind of river and watershed -improvement program 
requires some degree of alteration in the existing regimen of streams. 
The prevailing but generally erratic progression of floods and low- 
water flows may be changed by the building of impounding reservoirs, 
by diversions of water for beneficial uses, by soil-conservation measures 
and in other ways. The quantity of sediment transported by the 
stream may be changed as a result of deposition in reservoirs, by 
erosion-control measures on the watershed, or by revetment of the 
stream banks. The shape and slope of the stream may be changed 
by straightening, cut-offs, and jetties. Entirely new watercourses 
may be constructed to carry water diverted for irrigation, to provide 
drainage, or to create new navigable channels. 

If a stream is flowing in an alluvial valley over a bed composed 
mainly of unconsolidated sand or gravel, it is probable that the stream 
and its channel are essentially in equilibrium. The size and shape 
and slope are adjusted to the amount and variation of discharge and 
the supply of sediment of those sizes that make up its bed. If, then, 
some artificial change is made in the flow char act eristics, sediment 
supply or shape and slope of the channel, the stream will tend to make 
adjustments to achieve a new equilibrium. It will do so by scouring 
or filling its bed, widening or narrowing its channel, increasing or 
decreasing its slope. 

One of the most difficult problems encountered in open-channel 
hydraulics is the determination of the rate of movement of bed 
material by a stream. The movement of bed material is a complex 
function of flow duration, sediment supply, and channel character- 
istics. If a method is available for determining with reasonable 
accuracy the bed material movement under existing conditions, it 
would then be possible also to determine what the movement should 
be with any of these conditions altered. This would provide a 
reasonable basis for predicting what changes can be expected in a 
channel under a new set of conditions. 

Prediction of future channel changes has a very great economic 
importance in river-basin planning and development and in the 
operation and maintenance of river-basin projects. For example, if 
a large dam is constructed on an alluvial-bed river, all of the bed 
sediment normally transported will be trapped. The clear water 
released will tend to erode the channel bed downstream from the dam 
until a new equilibrium is established. Severe bed erosion may 
undermine costly installations such as bridge piers, diversion struc- 
tures, sewer outlets, and bank-protection works. 

Advance knowledge of the scour expected may influence the eleva- 
tion of tailwater outlets in power dams, influencing the power capacity 
of the dam to a very significant degree. On the other hand, regulation 
of the flow effected by the dam, particularly reduction in peak dis- 
charges, may make it impossible for the flow further downstream to 
transport all of the bed sediment delivered to the channel by tributary 
streams. Such a condition would result in aggradation of the main 




ity and adversely affecting 
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Differential reduction of peak flows and bed sediment supply by- 
watershed treatment measures — say 20-percent reduction in the former 
and 70-percent reduction in the latter — might initiate a cycle of 
damaging channel erosion in headwater tributaries and even down into 
main streams. If a channel is now aggrading, reduction in bed- 
sediment supply without proportional reduction in stream flow may 
be beneficial. If the channel is degrading, reduction in peak flows 
with less control of the sediment supply may be helpful. Except in 
areas where streams generally flow in rock-bound channels, the 
problems of bed-load movement and channel stability are almost 
universally present. Often they are critical if not deciding factors 
in not only the design and maintenance of works of improvement, but 
even of their feasibility. 

This publication does not attempt to give the specific solutions for 
all sediment problems in alluvial channels. It attempts only to 
provide a tool which the writer hopes is sufficiently general to apply 
to a large number of such problems. This tool is the bed-load func- 
tion. The equation for the bed-load function of an alluvial channel 
permits calculation of the rates of transport for various sediment sizes 
found in the bed of a channel which is in equilibrium. These equi- 
librium rates will be shown to be functions of the flow discharge.^ 

The significance of these equilibrium rates becomes apparent when 
one recognizes that they must have prevailed for a long time in order 
to develop the existing channels. By application of the bed-load 
function to an existing channel, it is possible to estimate the rate 
of bed-sediment supply. On the other hand, the same method may 
be used to determine the interdependent effects of changes of the 
channel shape, of the sediment supply, and of the flows in the channel. 

With the bed-sediment transportation rate a function of the dis- 
charge, it is clear that the long-term transport, that is, the average 
annual transport, can be predicted only if the long-term flow rates 
can be predicted. It will be shown that most sediment problems can 
be solved satisfactorily if at least the flow-duration curve is known. 
This fact emphasizes the urgent need for more knowledge about 
flow-duration curves for river sections of various sizes, for various 
climates, and for various watershed conditions. Today, accurate 
sediment-transport determinations are hampered more by a lack of 
necessary hydrologic data than by any other single factor. 

APPROACH TO THE PROBLEM 

The term "bed-load function" has proved to be useful in the descrip- 
tion of the sediment movement in stream channels. It is defined as 
follows: The bed-load function gives the rates at which flows of any 
magnitude in a given channel will transport the individual sediment 
sizes of which the channel bed is composed. 

This publication describes a method which may be used to deter- 
mine the bed-load function for many but not all types of stream 
channels. It is based on a large amount of experimental evidence, 
on the existing theory of turbulent flow, and beyond the limits of 
existing theory, on reasonable speculation. First, the physical char- 
acteristics of the sediment transportation process will be described. 
Next, sediment movement will be considered in the light of flume 
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experiments which allowed the determination of some universal con- 
stants of the various transportation equations. Finally, the calcula- 
tion of the bed-load function for a stream channel will be outlined to 
demonstrate the practical application of the method to determine 
rates of bed-load transportation. In its present state, the method 
appears to be basically correct. Although in various respects it is 
still incomplete, it appears to be useful for the solution of a considerable 
range of highly important problems. A special effort is made, 
however, to point out the unsolved phases of the problem. 

Some terms which recur frequently in this publication are defined 
as follows: 

Bed load: Bed particles moving in the bed layer. This motion 
occurs by rolling, sliding, and, sometimes, by jumping. 

Suspended load: Particles moving outside the bed layer. The 
weight of suspended particles is continuously supported by the fluid. 

Bed layer: A flow layer, 2 grain diameters thick, immediately 
above the bed. The thickness of the bed layer varies with the 
particle size. 

Bed material: The sediment mixture of which the moving bed is 
composed. 

Wash load : That part of the sediment load which consists of grain 
sizes finer than those of the bed. - 

Bed-material load: That part of the sediment load which consists 
of grain sizes represented in the bed. 

Bed -load function: The rates at which various discharges will trans- 
port the different grain sizes of the bed material in a given channel. 

Bed-load equation: The general relationsliip between bed-load rate, 
flow condition, and composition of the bed material. 

LIMITATION OF THE BED-LOAD FUNCTION 
The Undetermined Function 

Functions often become constant or even equal to zero under a wide 
range of conditions. That functions may not have any value in cer- 
tain ranges of conditions is mathematically demonstrated wherever 
the solution of the equation which defines the function becomes imagi- 
nary. But functions that become indeterminate under a wide range of 
conditions seem to be rather unusual. Unfortunately, the bed -Toad 
function has this character. 

In order to better understand this condition consider an example 
from the game of billiards. A player shoots the cue ball with the 
intention of hitting the red ball which is at rest. In what direction 
will the red ball move after the collision? Mathematically, the 
problem may be described in the following way: The independent 
variables are the angle a with which the cue ball rolls, and the original 
distance I between the two balls. The angle a, however, cannot be 

Predicted with mathematical accuracv, but is endowed with an error, 
'he actual value of a for any actual shot is defined by the angle oq 
which the player intends to use and by a small but absolutely random 
uncertainty a' which is only statistically determined by the accuracy of 
the player. The angle y at which the red ball begins to move after 
the impact depends upon the direction of the common plane of tangency 
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for the two balls at the moment of collision. Assuming no friction 
between the balls with a diameter D, the angle 7 may be calculated 
in terms of a', for instance, for an intended head-on collision if a and 7 
are measured from the original common centroid of the balls by the 
equation: 

sin 7 I 
sin a'~D 

As long as I is of the same order of magnitude as D, 7 will be of the 
same order as a' and it may be predicted with about the same ac- 
curacy as a. As soon as / becomes large compared to D, however, sin 
7 will be rather large and 7 may not be predicted with any degree of 

accuracy. With I larger than a given limit ^ may even become 

. larger than unity and the player may miss the ball completely. In 
this case the prediction of 7 from the intended average value <xq 
becomes meaningless because the possible fluctuations are much larger 
than this value itself. We may thus conclude that beyond a certain 
distance I the player, characterized bv an uncertainty a', has no 
chance at all of predicting 7 although he is able to do so with reasonable 
accuracy for small distances I. 

This example may show in a general way that physical problems 
exist which are determinate in part of the range of their parameters 
but are indeterminate in some other ranges. The transition between 
the two is usually gradual. The relationship between flow and sedi- 
ment transport in a stream channel is basically of this character. 
The critical parameter deciding the significance of the function in a 
given flow is the grain diameter of the sediment. 

The Alluvial Stream 

To introduce in simplified form the general case of sediment move- 
ment, assume a uniform, concrete-lined channel through which a 
constant discharge flows uniformly. Sediment is added to the flow 
at the channel entrance. Experience shows that sediment up to a 
certain particle size may be fed into such a flow at any rate up to a 
certain limit without causing any deposits in the channel. For all 
rates up to this limit, the channel is swept clean. An observer who 
examines the channel after the flow has passed can state only that 
the rate of sediment flow must have been below tins limiting rate; 
that is, below the "sediment transporting capacity" of the concrete 
channel. The sediment has not left any trace in the channel. Its 
rate of transport need not be related in any way to the flow rate. 
This kind of sediment load has been called "wash load" because it is 
just washed through the channel. 

If the rate of sediment supply is larger than the capacity of the 
channel to move it, the surplus sediment drops out and begins to cover 
the channel bottom. More and more sediment is dropped if the supply 
continues to exceed the capacity until the channel profile is sufficiently 
changed to reach an equilibrium whereby at every section the transport 
is just reduced to the capacity value. Now, an observer is able to 
predict that during the flow, sediment was transported at each section . 
at a rate equal to the capacity load; because if it had been more the 
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surplus would have settled out, and had it been less, the difference 
would have been scoured from the available deposit. Such a river 
section which possesses a sediment bed composed of the same type of 
sediment as that moving in the stream is called an "alluvial 
reach" (IS). 3 It is the main purpose of this paper to show how the 
capacity load in such an alluvial reach may be calculated. 

The Sediment Mixture 

This problem is highly complicated by the fact that the sediment 
entering any natural river reach is never uniform in size, shape, and 
specific gravity but represents always a rather complex mixture of 
different grain types. It has been found experimentally that the shape 
of the different sediment particles with few exceptions is much less 
important than the particle size. The specific gravity of the bulk of 
most sediments is also constant within narrow limits. It is, therefore, 
generally possible to describe the heterogeneity of the sediment mix- 
ture in a natural stream by its size analysis, at least when the mixture 
consists of particles predominantly in the sand sizes and coarser. As 
the derivations which follow do not introduce any molecular forces 
between sediment particles, they are automatically restricted to the 
larger particles, in general to those coarser than a 250-mesh sieve 
(Tyler scale) or 0.061 millimeters in diameter. This restriction does 
not seem to be serious, however, as most alluvial stream beds in the 
sense of the above definition do not contain an appreciable percentage 
of particles below 0.061 millimeters in diameter. 

Consider now how it may be necessary to modify the previous 
definition of the alluvial reach, of the sediment capacity, and of the 
bed-load function in view of the fact that all natural sediment supplies 
are very heterogeneous mixtures. Again, begin with the assumption 
of a flow in a concrete channel. Assume a sediment supply at the 
upper end of the channel, consisting of all different sizes from a maxi- 
mum size down through the silt and clay range. If the maximum 
grain is not too large to be moved by the flow, the channel will again 
stay clear at low rates of sediment supply. But an increase of the 
supply rate will eventually cause sediment deposition. 

Under most conditions only the coarse sizes of sediment will be 
deposited. It is true that a small percentage of the fine sediments 
may be found between the larger particles when the flow is past, but 
this amount is generally so small that one is tempted to conclude that 
these small particles are caught accidentally between the larger ones 
rather than primarily deposited by the flow itself. This is also sug- 
gested by the fact that the volume of entire deposit does not change 
if the fine particles are eliminated from it: they merely occupy the 
voids between the larger grains. 

A direct proof of the insignificance of these fine particles in the 
deposit can be found experimentally. The rate of deposit of the coai-se 
particles is a distinct function of the rate of supply. If more coarse 
sediment is supplied at the same flow, all this additional supply is 
settled out, leaving constant the rate which the flow transports through 
the channel. If only the rate of the fine particles is increased, how- 
ever, the rate of deposit of these particles is not influenced at all. 

8 Italic numbers in parentheses refer to Literature Cited, p. 68. 
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This basically different behavior of the fine and the coarse particles 
in the same channel has led the author and collaborators (<5) to assume 
that the fine particles in the flow still behave like material called "wash 
load" in the concrete channel, whereas the coarse particles act like the 
sediment in a strictly alluvial channel. These investigators give the 
limiting grain size between wash load and alluvial or bed load in terms 
of the composition of the sediment deposit in the bed. They state 
that all particle sizes that are not significantly represented in the 
deposit must be considered as wash load. More specifically, the 
limiting size may be arbitrarily chosen from the mechanical analysis 
of the deposit as that grain diameter of which 10 percent of the bed 
mixture is finer. This rule seems to be rather generally applicable as 
long as low-water and dead-water deposits are excluded from the bed 
sediment. 

Needless to say, the assumption of a sharp limit between bed load 
and wash load must be understood as a convenient simplification of a 
basically complex gradual transition. Virtually nothing is known 
about this transition today. This fact becomes apparent when the 
question is asked: what bed composition can be expected to result from 
a known sediment load in a known flow? No positive answer can be 
given to this question today. 

Another factor influencing the bed-load function is the shape of the 
channel cross section. If this section is not influenced either struc- 
turally or by vegetation it is only a function of the sediment and of the 
flow. We have today no clear concept of how to analyze these rela- 
tionships rationally even though we seem to have some workable rules 
for expressing the influence of the shape of the known cross section on 
the rate of transport. 

After this general discussion it becomes possible to define the pur- 
pose of this publication more specifically as the description of a method 
by which the capacity of a known alluvial channel to transport the 
different grain sizes of its alluvial bed at various flows may be 
determined. 

HYDRAULICS OF THE ALLUVIAL CHANNEL 

From the definition of the alluvial reach it was concluded that the 
transport of bed sediment in such a reach always equals its capacity 
to transport such sediment. It is easy to conclude from this that the 
flow is uniform or at least nearly so. The open-channel hydraulics of^_ 
nonuniform flow or the calculation of backwater curves is not par- 
ticularly important in this connection. Where such calculations are 
necessary for channels that are very actively aggrading or degrading, 
they are based on the use of the Bernoulli Equation as it is applied to 
channels with solid beds. 

- 

The Friction Formula 

The hydraulics of uniform flow include basically the description of 
the velocity distributions and of the frictional loss for turbulent flow. 
The writer has found that in describing sediment transport the velocity 
distribution in open-channel flow over a sediment bed is best described 
by the logarithmic formulas based on v. Karman's similarity theorem 
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with the constants as proposed by Keulegan (11). He gives the 
vertical velocity distribution as: 

^=5.50+5.75 log 10 ^=5.75 log 10 (9.05 ^) (1) 
for smooth boundaries and: 

•^=8.50+5.75 log 10 (*//k)=5.75 log 10 (30.2 yfk,) (2) 
u * 

for hydraulically rough boundaries. The transition between the two, 
including the rough and smooth conditions, may all be combined in 
the form: 

^=5.75 log 10 (30.2 ^)=5.75 log l0 (30.2 |) (3) 

whereby x is given in figure 1 as a function of k s /d. 
Herein are: 

Uy t he av era ge poi nt velocity at distance y from the bed 

u*= ^Talst^^SeHg the sheer velocity (4) 

st the density of the water 

S e the slope of the energy grade line 

R the hydraulic radius 

g the acceleration due to gravity 

y the distance from the bed 

v the kinematic viscosity of the water 

k t the roughness of the bed 

x a corrective parameter 

A= kjx the apparent roughness of the surface 

b the thickness of the laminar sublayer of a smooth wall: 

(5) 

The Friction Factor 

Next, a definition must be given for the roughness, k s , in the case 
of a sediment surface. For uniform sediment, k s equals the grain 
diameter as determined by sieving. Comparative flume experiments 
have shown that the representative grain diameter of a sediment 
mixture is given by that sieve size of which 65 percent of the mixture 
(by weight) is finer. 

A sediment bed in motion usually does not remain flat and regular 
but shows riffles or bars of various shapes and sizes. These irregu- 
larities have some effect on the roughness of the bed. Both flume 
measurements (6*) and river observations (7) have shown that this 
effect is rather considerable and cannot be neglected. An analysis of 
a large number of stream-gaging data in various rivers (7) has led to 
the following interpretation. 
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Resistance of the Bars 

The writer ($) has described a method by which the influence of 
side-wall friction on the results of bed -load experiments may be esti- 
mated. The assumption was made that the cross-sectional area may- 
be distributed among the various frictional boundaries in such a 
fashion that each unit will satisfy the same friction formula with the 
same coefficients that would apply if the entire cross section had the 
same characteristics. A similar approach can be used to describe the 
friction along an irregular sediment bed . It is assumed that on such 
a bed friction develops in two distinctly different wavs: (1) along the 
sediment grains of the surface as a rough wall with the representative 
grain diameter equal to k a ; and, in addition, (2) by separation of the 
flow from the surface at characteristic points of the ripples or bars. 
This separation causes wakes to develop on the lee side of the bars, 
characterized by rollers or permanent eddies of basically stagnant 
water such as those observed behind most submerged bodies of suffi- 
cient size. This flow pattern causes a pressure difference to develop 
between the front and rear sides of each bar so that part of the flow 
resistance is transmitted to the wall by this shape resistance, i. e. by 
normal pressures. 

Again we may be justified in dividing the cross-sectional area into 
two parts. One will contribute the shear which is transmitted to the 
boundary along the roughness of the grainy sand surface. The other 
part will contribute the shear transmitted to the wall in the form of 
normal pressures at the different sides of the bars. These may be 
designated A' and A" respectively. Both types of shear action are 
more or less evenly distributed over the entire bed surface and act, 
therefore, along the same perimeter. Two hydraulic radii may be 
defined as R f =A : f /p b and K"=A"\ Vh where naturally R'+R"=R b , 
the total hydraulic radius of the bed . 

This entire procedure of division may appear to be rather artificial 
since both actions occur along the same perimeter. The significance 
of this division becomes apparent, however, when one recalls that the 
transmission of shear to the boundary is accompanied by a transfor- 
mation of flow energy into energy of turbulence. This energy trans- 
formation caused by the rough wall occurs at the grains themselves. 
This newly created turbulence stays at least for a short time in the 
immediate vicinity of the grains and, as will be shown later, has a 
great effect on the bed-load motion. The part of the energy which 
corresponds to the shape resistance is transformed into turbulence at 
the interface between wake and free stream flow, or at a considerable 
distance away from the grains. This energy does not contribute to 
the bed-load motion of the particles, therefore, and may be largely 
neglected in the entire sediment picture. This may explain why the 
division of the shear into the two parts u, and u". is of first importance. 
We define: 

> (6) 

898379°— 51 2 
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From this it is understandable that the velocity distribution near a 
sediment grain in the bed surface is described by equations 1 to 3 
whereby u. assumes the value of u.. The average velocity in the 
vertical may be determined according to Keulegan as: 



^=3.25 + 5.75 log 10 (^)=5.75 log 10 (3.67 ^f) (7) 

for a hydraiilically smooth bed, and: 

4=6.25 + 5.75 log l0 (^)=5.75 log 10 (l2.27 (8) 

for a hydraulically rough bed. Again, the entire transition between 
the two cases inclusive of the extremes may be expressed by: 

4=5.75 log 10 (l2.27 ^p)=5.75 log 10 (l2.27 ^ (9) 

Where x is the same function of k,jb' as given in figure 4 (in pocket, 
inside back cover), and 

8'=±M-' (io) 

A corresponding expression uju'l may be calculated, and this 
expression must be expected to be a function of the ripple or bar pat- 
tern, basically corresponding to equation (9). The ripples and bars 
change considerably and consistently with different rates of sediment 
motion on the bed. We will see later that the sediment motion is a 
function of a flow function of the type: 



s/ R'S 9 



(11) 



wherein 5, and s f are the densities of the solids and of the fluid, 
respectively, Z> 35 the sieve size in the bed material of which 35 percent 
are ^iner, and R' and S e are as defined previously. The expression 
u/u+ may thus be expected to be a function of It was found that 
^such a relationship actually exists for natural, laterally unrestricted 
stream channels as given in figure 5 (in pocket, inside back cover). 
Any additional friction, such as from banks, vegetation, or other 
obstructions must, naturally, be considered separately. 

The Laminar Sublayer 

The presence of a laminar sublayer along a smooth boundary has 
already been mentioned . The thiciaiess 5 of this layer has been given 
as: 

(5) 
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in which v is the kinematic viscosity and the shear velocity along 
the boundary. 

Within the laminar sublayer the velocity increases proportionally 
to the distance from the wall: 

V 



(12) 



and at the edge of the layer where y=5, the velocity is: 



(13) 



From this point on out, the velocity follows the turbulent velocity dis- 
tribution of equation (1) which has the same value as equation (13) 
at y= 8 as is shown in equation (14): 



u«= [5.50+5.75 log 10 (11.6)] u*=11.6 



(14) 



The entire distribution is shown in figure 3. For an explanation of 
this distribution the reader is referred to any standard textbook of 

u y = uJ/5.75 log (30.2 -y) 




Fioubb 3. — Assumed velocity distribution near the laminar sublayer along a 

hydraulically smooth wall. 
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fluid dynamics where the viscous forces are shown to be strong enough 
in the laminar layer to prevent the occurrence of turbulence, whereas 
outside of the layer 5 the turbulent action is so strong that viscosity 
effects may be neglected. The transition between the two phases is 
naturally gradual and not as abrupt as shown in figure 3, but since 
nothing is known about the character of this transition the sudden 
change from one regime to the other may be used. 

The Transition Between Hydraulically Rough and Smooth Beds 

Along a rough wall the distribution of velocities is basically different. 
Einstein and El-Samni (8) have shown experimentally that the theo- 
retical boundary from which the distance y of equation (3) must be 
measured lies 0.2k g below the plane which connects the most prominent 
points of the roughness protrusions. The roughness k, is thereby 
given by the grain diameter of which about 65 percent by weight of a 
sediment mixture is finer. It is generally known, and may be seen 
from figure 4, that the wall acts hydraulically rough if fc,/5>5. The 
laminar sublayer as calculated, for a smooth wall has, then, a thickness 
of 5<Ar,/5. With the roughness protrusions 0.2^ high, this laminar 
sublayer would be strongly dissected by the protrusions if it existed 
at all. In reality it does not seem to exist in these flows because the 
shear is transmitted to the boundary differently. In order to under- 
stand this mechanism one must interpret the boundary as a sequence 
of bodies submerged in the fluid. For sufficiently large Reynolds 
numbers (k g /8) of the flow around the individual grains, separation of 
the flow will occur and a low-pressure wake of still water develops on 
the downstream side of the grains. The resultant of the normal pres- 
sures has a significant component in the direction of the flow which 
very soon (with increasing k s l8) becomes so large that all viscous shear 
can be neglected. This is the reason why friction formulas for flow 
along hydraulically rough walls (A: s /5>5) do not contain the Reynolds 
number in any form. 

The local velocity distribution at the rough wall follows according 
to El-Samni's measurements, equation (3), as close as 0.1 k s from the 
theoretical wall; i. e. even between the roughness protrusions. No 
measurements closer to the wall exist, but the assumption that the 
turbulent velocity distribution exists all the way down to zero velocity 
seems to be as valid as any other. This point is at a distance y 0 from 
the wall which may be determined from equation (3), such as: 

£=0=5.75 log,. (30.2 f) 

to be: 

k, 

y ° 30.2x 

which becomes: 

Vo 30.2 

for hydraulically rough walls. 



(15) 
(16) 
(17) 
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The Velocity Fluctuations 

All velocities introduced so far are time averages. The different 
types of sediment motion cannot be described by these time averages 
only. Movement in both suspension and along the bed can be ex- 
plained only if turbulence is introduced. Turbulence is an entirely 
random velocity fluctuation which is superimposed over the average 
flow and which can be described today only statistically. 

The turbulence velocity at any point of the flow has the following 
qualities: (1) It usually has three finite components, each of which has 
a zero time average. (2) The velocity fluctuations are random and 
follow in general the normal error law. Its intensity may be measured 
by the standard deviation of the instantaneous value. (3) Wherever 
shear is transmitted by the fluid , a certain correlation exists between 
the instantaneous velocity components in direction of the shear and 
in the direction in which the shear is transmitted. (4) As shown in 
recent measurements by Einstein and El-Samni (#), in the immediate 
proximity of a rough wall the statistical distribution of velocity inten- 
sities must be skewed since the pressure variations are following the 
normal error law there. 

The characteristics indicated in the four previous statements may 
need some explanation. The first statement is the easiest to under- 
stand. One may visualize turbulence as a complicated pattern of 
long eddies similar to the twisters and cyclones in the air, but smaller, 
more twisted and intricately interwoven so that they flow in many 
different directions. If water moves past a reference point with an 
average velocity, the eddy velocities assume all directions according 
to the various directions of the eddy axes. The only exception to 
this rule of three-dimensionality are the points very close to a solid 
boundary. There the velocities in direction normal to the boundary 
for obvious reasons are smaller than the two others. 

In statement 2 the expression "random" needs some explanation. 
If the intensity of the velocity component in one direction is recorded, 
the resulting curve resembles the surface of a very choppy sea, One 
may conceive of a periodic pattern like the waves. If one tries, how- 
ever, to find the amplitude and wave length of this curve it is apparent 
that both characteristics change constantly in an absolutely irregular 
pattern. It can only be concluded that the curve is continuous, that 
no discontinuities exist in the velocity itself.. The standard deviations 
of the velocity components have been measured, mainly in wind 
tunnels, by the use of hot-wire anemometers, where the standard 
deviation values may be determined directly. To the writer's knowl- 
edge, no measurements have been made sufficiently close to the wall 
to show the deviations mentioned in statement 4. 

Statement 3 is the basis of the so-called "exchange-theory" of 
turbulence which is today the most important tool for the study of 
quantitative effects of turbulence. Let us assume that a shear stress 
exists in the flow under consideration and that correspondingly a 
velocity gradient exists in the same direction. If, for instance, this 
shear is the consequence of the bed friction in a flow channel, the 
average velocities are essentially horizontal, but increase in magnitude 
with increasing distance from the bed. This increase is the velocity 
gradient previously mentioned. If an exchange of fluid masses is 
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visualized in this flow between two horizontal layers of different 
elevation, this exchange may be described by the flow velocities at a 
horizontal plane between the two. In this plane the amount of flow 
going up will equal the amount going down for reasons of continuity. 

All water particles going up have a positive instantaneous vertical 
velocity whereas the velocity in the downward direction is termed 
negative. All these fluid masses which move vertically through the 
horizontal plane have a horizontal velocity at the same time. Let us 
call the horizontal velocity fluctuation positive if the velocity is higher 
than the average at that elevation, negative if it is lower. The 
important point is that all water particles moving down through the 
plane from above originate from a region of higher average horizontal 
velocity. There exists a tendency, therefore, for the horizontal 
velocity fluctuation to be positive whenever the vertical velocity is 
negative. Similarly, the tendency is for the horizontal velocity 
fluctuation to be negative when the vertical component is positive. 
The correlation coefficient, which is the integral of the product of the 
two instantaneous velocities over a given time divided by the product 
of the standard deviation of the two, thus has a tendency to be negative. 
Its value gives a measure of the vertical movement of horizontal 
momentum through the plane, which represents a shear stress. 

This exchange motion transports not only mass and momentum 
through the reference plane, but also heat and dissolved and sus- 
pended matter, as explained under "Suspension." 

Statement 4, pertaining to the statistical distribution of static pres- 
sure near the bed, is based on empirical findings,, the significance of 
which is so far neither fully understood by itself nor in connection 
with the creation of turbulence in the boundary region. In this study 
the result has been used to great advantage, however, despite the lack 
of a full understanding of its general significance. 

SUSPENSION 

The finer particles of the sediment load of streams move predom- 
inantly as suspended load. Suspension as a mode of transport is 
opposite to what Bagnold called "surface creep" and to what he 
defines as the heavy concentration of motion immediately at the 
bed. In popular parlance this.has been called bed load, although as 
defined in this publication bed load includes only those grain sizes of 
the surface creep which occur in significant amounts in the bed. 

The characteristic definition of a suspended solid particle is that 
its weight is supported by the surrounding fluid during its entire 
motion. While being moved by the fluid, the solid particle, which is 
heavier than the fluid, tends to settle in the surrounding fluid. If the 
fluid flow has only horizontal velocities, it is impossible to explain 
how any sediment particle can be permanently suspended. Only if 
the irregular motion of the fluid particles, called turbulence, is intro- 
duced can one show that sediment may be permanently suspended. 

The effect of the turbulence velocities on the main flow was de- 
scribed in the discussion on hydraulics by reference to the fluid ex- 
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change. This same concept is used to describe suspension. Since the 
vertical settling of particles is counteracted by the flow, the vertical 
component of turbulence as described by the vertical fluid exchange is 
effective. Assume a turbulent open-channel flow. The section may be 
wide, the slope small. Consider a vertical sufficiently far from the 
banks to have two-dimensional flow conditions. On this vertical 
choose a horizontal reference section of unit area at a distance y 
from the bed. While the mean direction of flow is parallel to this 
area, the vertical velocity fluctuations cause fluid to move up and 
down through the section. Statistically, the same amount of fluid 
must flow through the area in both directions. To simplify the pic- 
ture, assume an upward flow of velocity (v) in half the area and a 
downward flow of the same velocity (— v) in the other half. The 

exchange discharge through the unit area is q v = ^v. If the exchange 

takes place over an average distance of l e at elevation y it can be 
assumed that the downward -moving fluid originates, as an average, 

from an elevation (y+^e) while the upward -moving fluid originates 

from (y— jj-iA The important assumption is made that the fluid 

preserves during its exchange the qualities of the fluid at the point of 
origin. Only after completion of the exchange travel over the dis- 
tance l e will jt mix with the surrounding fluid. From this it is possible 
to calculate the transport of a given size of suspended particles with 
a known settling velocity v s , if the concentration of these particles at 
y is Cy. The upward motion of particles per unit area and per unit 
time is: 

•HO"?*'" 0 

and the rate of downward motion is: 
The net upward motion is therefore: 

H^0 (l ' _, '' )_ 2 VIO 0 '^ (19) 

Neglecting all higher terms, the concentrations may be expressed as: 

1 , dc y 

C {.-\'.)- C "~2''JVj 



I l , (iCy 



(20) 
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Introducing these expressions (20) in equation (19) the net upward 
motion is: 

c * v ~2 vle dy 

Most interesting is the equilibrium status at which there exists no 
net flow in either direction: 

c ' v -+l v! ' § =0 (21) 

In equation (21) both v and l e are unknown. It is customary to 
assume that these two values are equal to the corresponding values 
in a similar equation for the exchange of momentum through the same 
area. Assuming that the shear due to viscosity may be neglected 
compared with that due to momentum transport, the depth d may be 
introduced: 

d—y_\ 

Ty— T 0 " "I — o 



•*[vm- i vm] 

[(%- \ ^y^+k 1 - §)] 



1 , du 
2 V,i 'Ty 

From this we calculate: 

2 Vt <~ Sf <* * d 'duldy 

Using equation (3) for the velocity distribution we may calculate 
dujdy: 



duy = 5.75 u± 
dy 2.303 y 

Introducing this value into equation (22): 



(23) 



i vl.=-0A0yuJ^2 (24) 



This value may be used in equation (21): 



c„v,=0A0yu, d ^ jjj (25) 
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Separating the variables: 

dc v= v, d dy . 
Cy 0.40u* y (d-y) K ' 0) 

and introducing the abbreviation: 

*=dk < 27 > 

we can integrate this equation from a to y 

r T* = / ^ a ° ge c ^ =l °^ ^)-log e (c a )=\og e (fy (28) 

' -f. a-r^ fe^ 

This may be rewritten in the form: 



and may be used to calculate the concentration of a given grain 
size with the settling velocity v, at the distance y from the bed, if the 
concentration c a of the same particles at distance a is known. 

The Transportation Rate of Suspended Load 



This entire derivation is based in part on the assumption that the 
instantaneous velocity of any suspended particle is that of the sur- 
rounding water plus its own settling velocity in this fluid, the two 
velocities being added vectorially. This makes the horizontal velocity 
component of the particle equal to that of the surrounding water. 
This allows us to calculate the flow rate of sediment particles at 
elevation y per unit area and time: c y u v . In order to see how this 
rate changes over the vertical we may assume Uy to be about constant 
as the logarithm changes very slowly compared to the power function 
of Cy, which equals zero at y=d and becomes infinite at y=o. An 
infinite concentration is an impossibility. None of the distributions 
can, therefore, follow equation (29) down to the bed, but there is no 
reason why they should not do so up to the surface. If the transport 
by suspension is integrated over the vertical it is very reasonable to 
begin the integration at the water surface and to integrate down to 
the depth y. We will see later how far down the suspension actually 
determines the transport. 

INTEGRATION OF THE SUSPENDED LOAD 

The integral of suspended load moving through the unit width of a 
cross section may be obtained by combining equations (29) and (3). 

£ c v u v dy=j d v c a j^jj 5.7 5u. log l0 (30.22//A)<fy (30) 

893379°— 51 3 
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This long expression may be slightly shortened by taking some of the 
constant factors out of the integral, and it may be simplified by 
referring the concentration to that at the lower limit of integration a. 
By replacing a by its dimensionless value A=a/d, and using d as the 
unit for y, we obtain 

2«=J c u u l/ dy=J dc v Uydy 

X'^Jlog.o^J (3.) 

In order to reduce the two integrals in equation (31) into a basic 
form, the log y is changed from base 10 to base e of the natural log- 
arithms using the relationship 

logic (y)=\og e (y) log 10 (e) (32) 

Aslog 10 (e) has the value of 0.43429 we may write equation (31) in the 
form 

L^ty+^LQ^)^***) (33) 

with 

v. 

2 = 



0.40 

y measured with d as unit 

A—aj d 

Herein are: 

q s the sediment load in suspension per unit of width, measured 
in weight, moving per unit of time between the water surface 
and the reference level y= a 
the shear velocity 

Ca the reference concentration at the level y=a. (c a is measured 
in weight per unit volume of mixture). 

A the dimensionless distance of this lower limit of integration 

from the bed. ^1=^ 
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z defined in equation (27) as the settling velocity v, of the 

particles divided by the Karman constant 0.40 and the 

shear velocity u* 
y the variable of integration, the dimensionless distance of any 

point in the vertical from the bed, measured in water 

depths d. 

Numerical Integration of Suspended Load 

Equations (33) and (34) are true to dimensions. Any consistent 
system of units will, therefore, give correct results. The two integrals 
in equation (33) cannot be integrated in closed form for most values of 
2. The numerical integration of the two integrals for a number of 
values of A and z was thus the only possible solution of the problem. 
After a survey of the available methods of approach it was decided 
to use the Simpson formula, integrating the two expressions in steps, 
whereby each series of integrations for a constant z value would pro- 
duce an entire curve of integral values in function of A. Each such 
integration was started from y=\ and proceeded toward smaller 
values of y. Table 1 gives a sample sheet for such a calculation. The 
values are calculated there for 2=0.6 and for 1>?/>0.1. By this 
same method the entire range l>i/>10" 6 was covered for the values 
2=0.2, 0.4, 0.6, 0.8, 1.0, 1.2, 1.5, and 2.0. In addition, 2=1.0 was 
integrated in closed form, as were some values for 2=0, 1.5, 2.0, 3.0, 
4.0, and 5.0. Values of z above 5.0, were considered unimportant for 
the problem in question, because only particles with very high settling 
velocities would have 2>5.0, and these particles move almost ex- 
clusively as surface creep. 

The calculations were then spot-checked by means of a method 
based on the development of the functions into binomial and poly- 
nomial series, with the original calculation carried through to 5 signi- 
ficant figures. It was found that the derived integral values were 
always correct to within at least 0.1 percent; i. e. to slide-rule ac- 
curacy. Most of the values given in table 2 are more accurate than 
can be obtained on a slide rule. 

Table 2 gives a list of all the values for the two integrals, J x and J 2 
calculated hy means of the Simpson formula. 

Table 3 gives in addition the comparison of some values as calculated 
by the Simpson rule with those determined by closed integration for 
the exponent Z=1.0. The other integrations checked in similar 
fashion with the largest deviation near A= 1 . for small values of z. 

Table 4 gives the integrals as solved in closed form for exponents 
Z=0, 3.0, 4.0, and 5.0. 

For practical use one needs many more values for the integrals than 
those calculated so far, even if the values of tables 2 and 4 cover the 
entire range of practically important A and z values. The full 
coverage of the entire field is accomplished more easily by graphic 
interpolation than by calculating additional integral values. For 
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this purpose it is convenient to transform equation (33) into the 
following form. 

fi.=n.6w {2.303 iog tt • Ji+ijj 

with 

y (35) 

Herein (11.6 u.) is the flow velocity at the outer edge of the laminar 
sublayer in case of hydraulically smooth bed, or the velocity in a 
distance of 3.68 roughness diameters from the wall in case of a rough 
wall. It is a good measure for the order of magnitude of the sediment 
velocitv near the bed. The symbol c a stands for the sediment con- 
centration at a distance a from the bed. The integral values Ii and 
I 2 are plotted in figures 1 and 2, respectively (see charts in pocket, 
inside back cover). 

LIMIT OF SUSPENSION 

Equation (29) shows that at the bed and in a layer near the bed 
the concentration becomes very high, infinite in the limit. Obviously, 
not more sediment particles can be present in any cubic foot than 
there is room for. The maximum is about 100 pounds per cubic foot. 
If, on the other hand, this concentration were assumed to exist 
at the wall, at y=0, then the rest of the section would have zero 
concentration. This further demonstrates that the suspended -load 
formula (29) cannot be applied at the bed. 

The mixing length l e has been defined as the distance which a water 
particle travels in a vertical direction before it mixes again with the 
surrounding fluid. The actual size of l e and of the fluid masses 
which as a unit comprise this exchange motion was not considered. 
It was pointed out, however, that l e is not a differential or a dimension 
of infinitely small magnitude. By correlation measurement, especially 
between velocities in wind-tunnel flows, it has been found that the 
order of magnitude of l e is the same as that of fluid masses which 
as a unit make the exchange motion. It was found, furthermore, 
that both decrease proportionally with the distance from the wall. 

The Bed Layer 

What happens to the suspended particle as it moves near the bed? 
Suspension, as defined herein, is obviously meant to describe the 
motion of a small particle moving around in a fluid with a rather small 
velocity gradient: Only then does the "velocity of the surrounding 
fluid" have a meaning. By means of the parameter l e this normal 
case of suspension may be easily expressed where the grain diameter 
D is very small compared to the mixing length l e . It has been shown 
that l e becomes smaller and smaller as the bed or any wall is ap- 
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proached. There, the velocity of the surrounding fluid as such has 
little meaning, because on one side of the grain the instantaneous 
velocity fluctuation may have one size and direction, and on the 
other side it may be entirely different. The resulting effect on the 
particle may be stresses in all directions, but the resulting force 
becomes negligible. One might say that the water only tickles/ 
the grain; it cannot push it. The value l e cannot be used there 
to express the particle motion according to equation (21). In the 
extreme case of very small l e values the particle is never supported 
by the flow and settles out. 

The flow layer at the bed in which the mixing length is so small 
that suspension becomes impossible has been found to be about 2 grain 
diameters thick. This may be designated as the "bed layer." In 
reality, the region in which suspension degenerates is not sharply 
defined. There exists rather a gradual transition to the rest of the 
flow. It is feasible, nevertheless, to idealize the condition by intro- 
ducing a sharp division between the bed layer and the bulk of the 
flow, as is customary in the case of the laminar sublayer. These two 
layers, however, are entirely independent, as indicated by the fact that 
the thickness of the bed layer is defined in terms of the grain diameter, 
which does not influence the sublayer. (See equation (5)). Since 
most sediment beds are mixtures of various grain sizes it is necessary 
to introduce a separate bed-layer thickness For each grain size, while 
the laminar sublayer is naturally the same for every grain in the mix-' 
ture. Certain difficulties arise when the two layers have about the \ 
same thickness. These difficulties will be covered later in discussing | 
the bed-load motion or surface creep which is the mode of motion 
inside the bed layer. 

PRACTICAL CALCULATION OF SUSPENDED LOAD 

The following equations are used for the calculation of suspended 
load: 

Equation (3) for velocity distributions. 
Equation (9) for the frictional loss along the bed surface. 
Figure 5 for additional shape resistance of sediment beds in natural 
rivers. 

Equation (29) for calculating the sediment concentration. 

Equation (34) for the total suspended sediment load. These equa- 
tions include a number of variables and parameters which need some 
explanation. First, equation (3) gives the velocity distribution near 
an even bed of the roughness k, (measured as a representative grain 
size). It may be applied, therefore, near the back of a bar. The 
value tt* assumes then the value v!^ The parameter z is to be read 
from figure 4. The roughness k s for a natural sediment mixture may 
be determined from the mechanical analysis of the bed material as 
that size of which about 65 percent (by weight) is finer. This rule of 
thumb has been checked very often and has never failed. 

Equation (9) and figure 5 need no additional explanation. 

In equation (29) the relative position of y and a is immaterial. 
Neither of them may be inside the bed layer, however, so that 2D is 
the minimum value for both.. The relation a=2D may be used as a 
rule of thumb. In the parameter z the settling velocity v, for sand 

898370°— 51 4 
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and gravel may be read from figure 6 (in pocket inside back cover), 
which represents a replot of Rubey's equation (17) and which seems 
to be fairly reliable, more so than theoretical curves based on the 
settling of spheres. Again, the value u* may be used for u* . 

It is clear from the preceding explanations that the different sedi- 
ment sizes must be calculated individually. In this connection the 
question arises as to how large the individual size ranges can be made 
without impairing the accuracy of the results. No systematic study 
has been made in this respect, but it was found that the convenient 

range of V2 as suggested by the standard sieve-sets is very satisfac- 
tory. Several times, when the range of sediment sizes in a river was 
rather large, the range according to sieves at a factor 2 was tried and 
no deviations larger than 10 percent were found. As a rule it will be 
sufficient to cover the bed sediment of a stream by six to eight size 
classes. 

Numerical Example 

All these relationships, graphs and rules may best be explained by 
an example. In a wide stream with a water depth of 15 feet and a 
slope of 2 feet to the mile a suspended load sample was taken 1 foot 
above the bed. The sample contained 1,000 parts per million of 
sediment of which 10 percent was coarser than sieve No. 60 (0.246 
millimeter) and finer than sieve No. 42 (0.351 millimeter). The 
representative roughness of -the bed is 1.0 millimeter and Z? 35 =0.5 
millimeter. With these figures in mind three questions may be 
asked : 

(1) How much sediment is moving in suspension above the sam- 
pling point per foot of width? 

(2) What is the concentration at the edge of the bed layer? 

(3) What is the total suspended load per foot of width? 
The solutions to these questions are found as follows : 

1. In solving the problem some preliminary values must first be 
calculated. 

If R' is assumed to equal R=d = 15 feet, the value of may be 
calculated using equation (11): 

T/ /2.68-l\ 0.5 5280 n AOA ■ * , , 

¥'=f — y~q — / 305 15^2 usmg feet, seconds, and 

pounds as the units. 

The following numerical values are used: £,=2.68; 1 mile= 
5,280 feet; and 1 foot =305 millimeters. 

u 

Figure 5 gives then -^=110 



Equation (8) gives u=y 15 32.2-5.75 log 10 ( 12.27 ^yjp)= 

11.7 ft. /sec. with the acceleration of gravity assumed to be 
<7=32.2 ft./sec. 2 
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Now u* may be calculated u% = -~q = 0 . 1 0 6 foot per second . From 

*u nr/ m 1 0.106 2 5 2 80 A AOA . 
this ^= 2.32.2 =0 - 920 f00t ' 

This gives a first correction for R' 

#'=15.0-R"=14.08feet. 

Now a second approximation is obtained: 

2,68-1 OS 5280 _J_ 
V ~ 1 305 2 14.1~ U0A0 

4r=98 



^=V 141 5^0 32,2 5 * 75 lQ gio(l2.27 ^ 305) 
= 11.3 feet per second 
113 

14'=-^- =0.1 15 feet per second 

Straight line extrapolation permits us to try: 

£"=1.13 feet 
5'= 13.87 feet 



^V 13 ' 87 5^0 32,2 5.75 log 10 (l2.27i^ 305) 

= 11.2 feet per second 

112 

ul'=-x^- = 0.1 19 feet per second 

R" = ( **oo = 1.16 feet, which is sufficiently close 

#' = 15.0-1.16=13.84 feet 



C = ^/l3. 



84 rr^r 32.2 = 0.411 feet per second. 
5 2oU 
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The average grain size D for the range is given by the geometric 
mean of the range limits 



D= V0.351 . 0.246=0.294 mm.=0. 000964 feet. 
At the point of sampling ^1=^^=0.0667 

k,= \ millimeter =0.00328 feet 

kg k $ u* = 0. 00328 0.411 
d~ll.§v~ 11.6 0.000012 

z=1.00 (from figure 1) 

C=0.10 J-^y^oo 62.4 = 0.00624 pounds per cubic foot. 

The exponent z is calculated from the settling velocity 

a, =0.1 25 feet per second 

v, _ 0.125 
2 ~0.40 <~0.40.0.411~ U/OU 

using equation (35) and the graphs for I x and I 2 , figures 1 and 2, 
respectively 

7,=0.59 

-0.87 

q,= 11.6 ■ 0.411 • 0.00624 • ^{^^(o^^ 15)(0.59)-0.87 j 

=0.0298(6.99-0.87} =0.183 pounds per second-foot 

which is the sediment transport per foot of width above the 

sampling station. 
2. The sediment concentration at the ed^e of the laminar sublayer 
for sediment between 0.246 and 0.351 millimeters is calculated from 
equation (29) directly. 

y =2Z>=2-0.000964=0.001928 feet 

a =1.0 foot 

Cb— 0.00624 pounds per cubic foot 
z =0.760 ✓ 
d =15.0 feet . 
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m fd-y a V ftftnft9A / 15.0 1.0 V ™ 
C ' =C ' \~J- d=i) =0 00624 V0^01928 Tit) 

=0.00624 (555) 0 - 760 

=0.764 pounds per cubic foot 

3. The total suspended load of particles between 0.246 and 0.351 
millimeters is calculated from equation (34) for the reference level 
A at a=2D. 

10 

c a = 0.764 pounds per cubic foot 
2 =0.760 

*ig 7 j* 1 " 0111 g^P* 18 * figures 4 and 5, respectively 
g a =11.6.0.411.0.764.0.001928|^log 10 (^|^_15)56.2-19.7| 

=0.00702[66.5— 19.73=0.329 pounds per second-foot of which 
0.329—0.183=0.146 pounds per second-foot move within a foot 
from the bottom. 

This provides the full solution of the problem. The only part which 
cannot be calculated directly is the division of R into R' and R". 
Trial-and-error methods must be used. It is possible, however, to 
calculate the entire curve at R against R' by first assuming values of 
R', then calculating R" and 7? and plotting the results as a curve. 
The value of R' may then be read from the curve for any R value. 

BED-LOAD CONCEPT 

It has been demonstrated that the motion of sediment particles in 
the bed layer cannot be described by the theory of suspension. The 
reason is that the particles there are not "suspended" by the fluid. 
They settle out, down to the bed. This does not imply, however, 
that they do not move any more. It only means that, while moving, 
their weight is supported by the noninoving bed and not by the fluid. 
Accordingly, they move by rolling and sliding on the bed or by making 
short hops (of a few grain diameters in distance), more or less con- 
tinuously remaining in the bed layer while moving as bedload or 
surface creep. The expression "surface creep" is more frequently 
used for particles that move at times in suspension. For the large 
particles which never go into suspension in the flow the more familiar 
expression "bed load" may be used without any danger of misunder- 
standing. 

Bed-load motion has been studied principally in laboratory flumes 
under conditions where suspension may be neglected. An exception 
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is the classic series of experiments by Gilbert (9) which included sizes 
down to 0.3 of a millimeter. A large number of authors have tried 
to find a single bed-load equation which would describe the total 
transport in all these experiments. The writer believes that this is 
basically impossible as movement in suspension follows principles 
which are entirely different from those which govern bed-load trans- 
port. Even if a general formula could be devised to describe all of 
Gilbert's experiments it is almost certain that this formula would not 
be applicable to the great water depths of a natural river. 

For the present purpose, bed load moving as surface creep is defined 
strictly as: the motion of bed particles in the bed layer by rolling, 
sliding, or hopping. This definition purposely excludes from bed load 
all particles finer than those of the bed. For practical purposes one 
may even exclude the finest 10 percent (by weight) of the bed material 
since these particles do not usually represent a structural part of the 
bed but only loosely fill the pores between the larger particles. 

The qualitative results of two flume experiments which can be 
easily repeated may be helpful in understanding the mechanics of 
bed-load motion. In the first experiment a flow is discharged con- 
tinuously over a sediment bed and sediment is added at the upper 
end until deposition causes equilibrium to be established throughout 
the length of the bed. Then certain particles, marked so they can 
be identified, are added at the upper end. 

Visual observation shows that the bed-load particles move with 
a velocity that is comparable with the water velocity near the bed 
made visible by injecting dye. By assuming that the bed particles 
move at the same velocity as the flow in the bed layer, the time after 
which the marked particles should have reached the downstream end 
of the flume can be calculated. If 100 percent is added to this time 
for a safety allowance, one might expect to find all the marked sedi- 
ment particles safely in the deposit at the downstream end of the flume. 
On the contrary, however, if the flow is interrupted at that instant 
and the deposit inspected one may find only one or two of the marked 
particles there. Most of them have traveled only a small fraction 
of the distance and are found in the stream bed near the upper end 
of the flume. This result is not compatible with the assumption of 
an equilibrium condition unless an equal number of bed particles 
have been scoured from the bed during the same period. This possi- 
bility may be tested by a second experiment. Before the experiment 
is begun the water is drained from the flume and the bed allowed to 
dry. Dye is sprayed on a predetermined part of the bed area and 
thus all the sediment particles of that area are marked. Upon resump- 
tion of the experiment one observes that, gradually, all marked particles 
are eroded and replaced by others of the same type. 

The two experiments prove definitely that an intimate relationship 
exists between the bed-load motion and the bed. Actually, bed-load 
motion is motion of the bed particles. In experiments similar to the 
one described first, the author (#) has shown that the motion of the 
bed particles is fully governed by statistical laws which can be stated 
as follows: 

1 . The probability of a given sediment particle being moved by 
the flow from the bed surface depends on the particle's size, shape, 



Digitized by Google 



THE BED-LOAD FUNCTION FOR SEDIMENT TRANSPORTATION 31 

and weight and on the flow pattern near the bed but not on its 
previous history. 

2. The particle moves if the instantaneous hydrodynamic lift 
force overcomes the particle weight. 

3. Once in motion, the probability of the particle's being rede- 
posited is equal in all points of the bed where the local flow would 
not immediately remove the particle again. 

4. The average distance traveled by any bed-load particle be- 
tween consecutive points of deposition in the bed is a constant 
for any particle and is independent of the flow condition, the 
rate of transport, and the bed composition. For the sediment 
grain of average sphericity this distance may be assumed to be 
100 grain-diameters. 

5. The motion of bed particles by saltation as described by 
Bagnold (1) may be neglected in water, as proved by Kalinske 
(JO). 

6. The disturbance of the bed surface by moving sediment- 
particles (1) may be neglected in water. 

From these findings it follows that the variables which at any spot 
in the bed determine the bed load are (1) the composition of the bed 
within an area 100 grain-diameters from the spot, and (2) the flow 
conditions near the bed in this same area. Conversely, the bed-load 
rate is not influenced by a slow change of the bed location as long as 
the composition of the bed does not change. The laws of the equi- 
librium transport may be used, therefore, to describe the bed-load 
transport on a changing bed, as long as it is possible to describe the 
bed and the flow locally during the transition. 

Some Constants Entering the Laws of Bed-Load Motion 

The following derivation of a bed-load equation is based on a large 
amount of data obtained from flume studies and from field measure- 
ments and observations. One set of experiments determining the 
forces acting on bed particles (8) contributed the following informa- 
tion : 

1. The roughness diameter k, is somewhat larger than the 
average grain diameter. It may be assumed to equal the diam- 
eter Uu which is the sieve diameter of the grain of which 65 
percent of the mixture (by weight) is finer. 

2. Friction along a sediment bed without ripples, bars or other 
irregularities is well described by equations (3) and (8). 

3. Ad average dynamic lift force acting on the surface particles 
of the bed may be expressed as an average lift pressure p L by the 
equation ^ 

PL=c L s f ^ (36) 

where d,=0.178, s f the density of the fluid, u the flow velocity at 
a distance 0.35 Z) 35 from the theoretical b.ed, and Di& is the sieve 
size of the grains of which 35 percent are finer. 

4. The pressure fluctuations due to turbulence follow in their 
duration the normal error law, the standard deviation being 0.364 
of the average lift. 
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Besides these experiments there are the results from a large number 
of bed-load experiments which determined the sediment transport for 
a given flow over a given bed. 

The least complicated case of bed-load movement occurs when a 
bed consists only of uniform sediment. Here, the transport is fully 
defined by a rate. Whenever the bed consists of a mixture the trans- 
port must be given by a rate and a mechanical analysis or by an 
entire curve of transport against sediment size. For many years this 
fact was neglected and the assumption was made that the mechanical 
analysis of transport is identical with that of the bed. This assump- 
tion was based on observation of cases where actually the entire bed 
mixture moved as a unit. With a larger range of grain diameters in 
the bed, however, and especially when part of the material composing 
the bed is of a size that goes into suspension, this assumption becomes 
untenable. Some examples of such a transport are given in the flume 
experiments described on pages 42 to 44 of this publication. 

The mechanical analysis of the material in transport is basically 
different from that of the bed in these experiments. This variation 
of the mechanical analysis will be described by simply expressing in 
mathematical form the fact that the motion of a bed particle depends * 
only on the flow and its own ability to move, and not on the motion 
of any other particles. 

THE BED-LOAD EQUATION 

The bed-load equation by definition is the equation which relates 
the motion of bed material per unit width of bed layer to the local 
flow. After the description of motion in the preceding chapter it may 
be easily understood that this equation must express the equilibrium 
condition of the exchange of bed particles between the bed layer and 
the bed. For each unit of time and of bed area the same number of 
a given type and size of particles must be deposited in the bed as are 
scoured from it. 

To express the rate at which a given size of sediment particles is 
deposited in the unit bed area per unit of time, let q B equal the rate 
at which bed load moves through the unit width of cross section and 
let i B equal the fraction of q B in a given grain size or size range. Thus 
q B i B is the rate at which the given size moves through the unit width 
per unit of time. All the particles with a particular diameter D are 
just performing an individual step of 100 D or, more generally, of 
AjJ) length. When they pass through the particular cross section 
where q B is measured, however, it is not known what part of AjJ) the 
individual particles have already travelled. They must be assumed 
to be deposited anywhere from zero to A J) downstream of the section. 
The area of deposition is A J) long and has unit width. If q B is 
measured in dry weight per unit time and width and if A 2 D 3 is the 
volume of a particle, a, its density, and g the acceleration of gravity, 
the number of such particles deposited per unit time in the unit of 
bed area may be expressed as: 

* 

AlDA^^q A 2 A L gs s IP 
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The rate at which sediment particles of this size are eroded from the 
bed per unit of time is proportional to the number of particles exposed 
at the bed surface per unit of area and to the probability p, of such a 
particle being eroded during a second. If i b is the fraction of the bed 
sediment in the given size range it may be assumed that this repre- 
sents also the fraction of the surface covered by particles in the same 
size. The number of particles D in a unit area of bed surface is thus: 

AD 2 

and the number of particles eroded per unit area and time is: 

AD 2 

If the time, t u necessary to replace a bed particle by a similar one 
were known, the probability of removal p s per second could be re- 
placed by the absolute probability p to be exchanged as p t h=p. Thus 
it follows that p s is the number of exchanges per second, t x the time 
consumed by each exchange, and p s ti the total exchange time per 
second, or the fraction of the total time during which an exchange 
occurs, which is the definition of p. The number of particles eroded 
per unit area and time is then: 

ibp 
A X D% 

The Exchange Time 

No method exists today of determining experimentally the exchange 
time ti. But experiments indicate that t x is another characteristic 
constant of the particle like the unit-distance of travel. As such it 
must be possible to describe it without introduction of the flow. The 
time ti may then be assumed to be proportional to the time necessary 
for it to settle in the fluid through a distance equal to its own size: 

and the number of particles eroded per unit of area and time is: 

Up lobs—st) 

A 1 D 2 A 3 \ Ds t 

and the bed-load equation (38) shows that this rate of scour equals 
the corresponding rate of deposit: 

Hq.B _ UP jghs — *t) 

i^A^^A.AD 2 ^ Ds t ^ 

M0»»79° -51 5 
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The Exchange Probability 

The probability, p, of being eroded has been defined as the fraction 
of the total time during which at any one spot the local flow conditions 
cause a sufficiently large lift on the particle to remove it. With all 
points of the bed statistically equivalent, p may also be interpreted 
as the fraction of the bed on which at any time the lift on a particle 
of a ^iven diameter D is sufficient to cause motion. 

With this interpretation, p may be used to calculate the distance 
AJ) that a particle travels between consecutive places of rest. It 
has been shown that this distance AJ) was found empirically to be a 
constant for each size particle. As long as p is small, deposition of the 
particle is practically everywhere possible and A L equals a general 
constant, X, which has about the value 100. If p is not small, however, 
it must be recognized that deposition cannot occur on that part 
(p) of the bed where the lift force exceeds the particle weight. By 
averaging the distances traveled by the individual particles until they 
are able to settle out, the value AJ) can be expressed as: 
(1— p) particles are deposited after traveling \D 
p particles are not deposited after traveling \D. Of these, 
p (1 — p) particles are deposited after traveling 2\D 
p 2 particles are not deposited after traveling 2A/X Of these, 
p 2 (1 — p) particles are deposited after traveling 3XZ) and so on. 
The total (and average) distance traveled by the unit is obtained by 
addition : 

A L D = (1 -P)p"(n+ l)X/?= ( /f^ (39) 

as may be found easily. If this value is introduced in the above bed- 
load equation it may be rewritten: 

gfljfi(l-/>)^ ibP 1 9 '(us/) , 4m 
~A 2 \gs t D* AiA^D 2 \ Ds f 

or, separating p on one side of the equation : 

p VAiAfl /ibMI* ( * V ( 1 y \ 
1-p L ^ J \iJhs9 \u-tfj Kg&J J 

j^ ir [A,) l fW = A^ (41) 

Therefore, p is the probability of a particle being eroded from the bed 
and * is defined as: 

Thus * is a dimensionless measure of the bed-load transport; it 
may be called the intensity of bed-load transport. Being a dimension- 
less parameter it does not change with the scale and is, therefore, 



Digitized by Google 



THE BED-LOAD FUNCTION FOR SEDIMENT TRANSPORTATION 35 

invariant between model and prototype. This relation may also be 
expressed as follows: If 3> is equal in two different flows, the two rates 
of bed-load transport are dynamically similar. 

Determination of the Probability p 

As already noted, p is the probability of a particle being eroded 
from the bed, which means that the probability of the dynamic lift L 
on the particle is larger than its weight (under water). The weight of 
the particle under water is 

W'^gk-^Atl? (43) 
while the lift force may be expressed as 

L=c L s t \ v?A x B> (44) 

In these two expressions all variables have been defined previously 
except the lift coefficient c L which Einstein and El Samni (8) found by 
measurement to be Ci,=0.178, and the velocity u near the bed which 
El Samni found must be measured at a distance of O.'SoD from the 
theoretical bed for uniform sediment. 

The forces acting on individual particles of a natural sediment mix- 
ture in a bed cannot very well be measured. They must be deter- 
mined from their effect on the movement of particles. In analyzing 
the experiments described on pages 42 to 44, the following general 
results were found: 

L The velocity acting on all particles of a mixture must be meas- 
ured at a distance 0.35A from the theoretical bed , whereby : 

A'=0.77A if A/5>1.80 > | 

(45) 

X= 1.395 if A/6<1.80J 

2. The particles smaller than X (X^>D) seem to hide between the 
other particles or in the laminar sublayer, respectively, and their lift 
must thus be corrected by division with a parameter £ which itself is 
a function of D/X (fig. 7, in pocket, inside back cover). 

3. An additional correction factor Y was found to describe the 
change of the lift coefficient in mixtures with various roughness con- 
ditions. Figure 8 (in pocket, inside back cover) gives the correction 
Y in terms of k,/8, Y being unity for uniform sediment. 

Using these assumptions the velocity in the expression for the 
average lift L may be written as: 

u=^5.75 log, 0 ( 1 

u*=KS.g 5.75 2 log? 0 (l0.6 ^) 
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At any instant the lift force may be described by: 
' X=0.178s 1 ^ 1 Z? 2 ii?;^5.75 2 log? 0 (^10.6^)[l + 7 ? ] (46) 

where rj is a parameter varying with time. 

Now p may be expressed as the probability of W'/L to be smaller 
than unity: 

I W*r 1 1 U-H D ) ( 2A 2 \ 1 
' L Ll+r/Jl st V.0.178 .A, 5.75V log? 0 (10.6 X/A) 

(47) 

The value of ij in this inequality may be either positive or» negative. 
In both cases the lift is actually positive and must, therefore, be 
understood on an absolute basis. The inequality may be written 
in absolute values : 



|[1 + ,]|>U»^ 

Introducing the abbreviations: 
,-ft D 



(48) 



B- 



S t R' B S e 

0.178 A x 5.75 2 
^v=log 10 (10.6A7A) J 



(49) 



Introducing the two correction factors £ and Y according to the pre- 
vious!}' quoted assumptions, the inequality (48) may be generalized: 



l+il\>tY**Z- 

Px 

where 

£ is a function of D/X (figure 7) 1 
Y is a function of k t /8 (figure 8) 
B'=BJ0 2 
0=log io (10.6) 
and 

p 2 /p 2 x = l for uniform grain and x=l 
Y=l for uniform grain and x=l 
1=1 for uniform grain and x=l 



(50) 



(51) 



(52) 
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Inequality (48) may be written more conveniently by squaring and 
division by rj 0l the standard deviation of tj. Introducing 17 =170*?* 

[1/%+iJ 1 > Y*Bl*> (0*/&) a =Bl*i (53) 

and^ = $Y(/9 2 //3 2 A)^ ) • 
Using these symbols, the limiting cases of motion may be written as: 

or 

[ljjLIMIT = ±B+*+— If no (55) 

As the probability for 17* values is distributed according to normal 
error law, the probability p for motion is: 

p=\ — U\ e~**dt (56) 

where t is only a variable of integration. 

By combination with equation (41) the final bed-load equation is 
obtained: 

This equation appears to be very complicated and difficult to use, but 
it is rather easy to apply. ifo, A* and B* are universal constants such 
that the equation may be represented by a single curve between the 
flow intensity ^ + and the intensity of bed-load transport This 
relationship may be calculated from tables of the probability integral 
for the value 1/tjo=2.0 as determined by El Samni. The constants A+ 
and B+ were determined from bed-load experiments with uniform 
grain for which $*=* and = Figure 9, appendix, shows a 
plot of some experimental points with the curve using ,4* =27.0 and 
B* =0.156. As experiments two series were used: the low intensities 
are flume studies made by the writer in 1932-35 at Zurich, Switzerland, 
using gravel of about 27 millimeters average size in a 7X7 foot channel; 
the higher intensities are the Gilbert experiments with 0.785 millimeters 
average grain size. No measurements were left out from either set. 

Field experience with the applicability of these formulas is still 
limited. The formula can unquestionably be applied to coarse sedi- 
ment as it is almost identical with most other bed-load formulas for 
low intensities. For the higher intensities which occur only with small 
particle sizes, some applications to actual rivers have given encourag- 
ing results, while more applications under a wider range of conditions 
are still necessary to prove its universal applicability. No failure has 
been encountered to date, however. 
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Although the correction Y as a function of kj8 does not require any 
explanation, the correction factor |, which gives the effect on the 
transport when the small particles of the bed hide either behind and 
between larger particles or in the laminar sublayer, needs some 
comment. In calculating transportation rates for particles affected 
by £, one finds that the rate calculated on the basis of the curve of 
figure 7 rapidly becomes negligible as higher values of £ are approached ; 
i. e. as the grain becomes small compared to X. This result seems to 
contradict the general observations in river-sediment measurements 
to the effect that these fine particles actually represent the bulk of 
the total load. But this only seems to be a contradiction. In the 
earlier part of this publication, it was explained that the entire sedi- 
ment load of a stream can be divided into two parts: namely the bed- 
material load, for which the bed-load function may be established, 
and the wash load for which no such relationship exists. 

The division between the two parts of the load was made rather 
arbitrarily at a special grain size determined from the grain composi- 
tion of the bed. In reality, however, this sharp division between the 
two does not exist and there is instead a gradual transition. As a 
result, in this range of sizes part of the load is bed-material load, part 
is wash load. Or put another way, the rate at which such a border 
particle size moves cannot be less than the bed-load rate if the bed is 
not to be changed. Additional wash load may or may not move 
without any appreciable effect on the bed, and, therefore, the maximum 
possible loed which may move without causing a change of the bed by 
deposition can be considerably higher. 

In this connection one question of interest is why most flume 
experiments have given the nimimum transport without any wash 
load. The reason apparently must be found in the experimental 
procedure applied in the experiments. The bed was first filled into 
the flume and the load or transport represented eroded bed particles. 
According to the preceding explanation, that is merely the minimum 
load condition. An alluvial nver, on the other hand, often operates 
under entirely different conditions. Its bed is being maintained by 
deposition which equals or exceeds the scour. Under this condition, 
wash load may be expected to occur. Actually this condition has not 
so far been tested in flume studies, but there* is no reason why that 
could not be done. 

TRANSITION BETWEEN BED LOAD AND SUSPENDED LOAD 

Up to this point, suspension has been described by the exchange 
theory between different layers of a turbulent flow. The sediment 
load is described in the resulting formulas in the form of concentra- 
tions, or, to be more exact, by ratios of concentrations. Special 
emphasis was put on the fact that the suspended load theory permits 
determination only of the local distribution of the sediment but not 
of the absolute amount or rate of transport. Experience in river 
measurements, on the other hand, shows that certain parts of the 
suspended load follow a function oi the flow and represent, therefore, 
part of the bed-load function, as previously defined. That fact 
cannot be explained by the theory of suspension. Another relation- 
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ship must determine the sediment concentration at some reference 
elevation in the vertical. 

In a search for this relationship it is helpful to remember that it 
exists only for particle sizes which are represented in the bed of the 
stream (5). This fact suggests that the concentrations must be gov- 
erned from the bed up in some fashion. It has been shown how the 
flow together with the sediment composition of the bed determines the 
transport of sediment in the bed layer. The relationship was found 
to be governed by the exchange of sediment particles between the 
bed layer and the bed. It is, therefore, very probable that the rela- 
tionship governing the concentration at the lower edge of the suspen- 
sion could be found by setting up an expression for the exchange of 
sediment particles between the suspension and the bed layer. 

Such an equilibrium condition must have exactly the same char- 
acter as normal suspension. Actually, the normal suspended-load 
calculation can be extended down to the bed layer. Only the con- 
centration at the upper boundary of the bed layer must be determined. 
The condition of exchange can then be expressed by equating the con- 
centration at the upper boundary of the bed layer with the bottom 
concentration of the suspension above. 

The problem is now reduced to the determination of the concen- 
tration at the upper boundary of the bed layer. No available experi- 
mental data either support or contradict directly any assumptions in 
this phase of the sediment problem. The total rate of transport, 
(/ B i B , of a given grain size in the bed and the thickness of the layer 
of 2 D within which this transport occurs have been determined and 
assumed, respectively. From these values the average concentration 
in the bed layer may be found. The concentration is defined as the 
weight of solids per unit volume of water-sediment mixture. First, 
the weight of bed-load material in motion may be calculated for the 
unit of bed area. Let u B be the average velocity with which bed-load 
material moves in the bed layer while in motion, not including the 
rest periods. Then the weight of particles of a given size per unit 
area is: 



u B 

The volume of the unit area of bed layer is 2 D and the average con- 
centration in the layer is: 

u B 2.D 

It is probably not much in error to assume that the concentration in 
the entire bed layer is constant, since the layer is only two diameters 
thick. In order to leave open a possibility for correction, the following 
equation may be set up : 

"°-~2DW (58) 
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The velocity u B is not known. Both the flow velocity and the trans- 
port near the bed are functions of u*. The two must determine u B . 
This makes it very probable that u B is proportional to ui because it 
has the dimension of a velocity, too. Therefore, c a may be expressed 
in the form : 

*-*sifc < 59 > 

The value A 6 must be determined experimentally. It includes 
both the distribution of concentrations in the bed layer and the 
velocity of the bed load. It can, therefore, best be determined from 
flume experiments which compare the suspended load with the hydrau- 
lics of the flow. The set of experiments described further on under 

1 'Flume Tests With Sediment Mixtures" suggests that A i = T ^-z 
is an average value. Equation (59) may thus be written as: 

2s=11.6 c a v.* a (60) 

and the total suspended load per unit width, i t q t , may be calculated 
from equations (34) and (60) 

^• = ^|o^ loglo (l?) /i+Iij^WPZi+JJ (61) 
whereby 

has the same value for all different grain sizes of a section. This 
relationship relates transportation as bed load to that in suspension 
of all particle sizes for which a bed-load function exists. 
The total load q r may now be calculated 

i r g r =r B g B +;,g 5 =; / ,g B {P/i + / 2 +l} (63) 

in which, as shown earlier, I 2 is always negative. This completes the 
presentation of the theory on the basis of which the bed-load function 
of a reach may be calculated. 

THE NECESSARY GRAPHS 

The following comments further explain the use of the working 
graphs (figs. 1, 2, and 4 to 9, inclusive) referred to in the course of the 
preceding explanation of the theory, and figure 10, work graph, 
referred to below. 

1. Figures 1 and 2 give the two integrals I x and I 2 in terms of the 
exponent z and the limit A. Of these, I x is always positive and I 2 is 
always negative. For values of s>5 the expression P/i-f 1 2 usually 
becomes smaller than 0.2 and may be estimated directly as it may be 
practically neglected against the additional 1 in equation (63). This 
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means that the particles do not go into suspension but stay perma- 
nently in the bed layer. 

2. The curve of x against k 9 (6 of figure 4 is derived from Nikuradse's 
experiments (15, 16) , which used sand grains glued to steel pipes as 
roughness. The curve has always been found to describe reliably 
the roughness of plain sand beds. Some deviations may be expected 
if the range of grain sizes in the deposit is very large. Especially, 
the rather crude rule of k,= D K may not be too reliable if the sediment 
contains appreciable amounts of silt and clay. Z?m is the grain diam- 
eter of which 65 percent by weight is finer. The curve cannot be 
used for most fabricated surfaces such as steel, galvanized pipe, and 
pitted concrete. 

3. The curve of u/tu against & is less reliable than the x curve. 
It attempts to describe from existing river measurements the effect 
of the irregularities of a natural stream channel on channel roughness. 
Considering the wide variety of natural channels it is readily recog- 
nized that no accurate appraisal of this effect can be given. The 
curve given in figure 5 seems, however, to describe rather closely the 
behavior of natural channels not constricted by artificial banks, 
vegetation, or other obstructions. Channels dissected by trees, by 
stable vegetated islands, or by rock islands may show a u/u* value 
which is only about 0.7 of the curve value. Some reaches with ex- 
tremely heavy tree growth have showed only 0.5 of the curve value. 
The curve does not describe friction conditions in flumes. Flume 
measurements may come down to the curve, especially at high rates 
of transport, but they always seem to be much higher at low rates of 
transport. The side walls seem to straighten the flow so much that 
they prevent much of the channel irregularity from occurring. This 
is one of the most important reasons for the uncertainty of the curve: 
it can be developed only from river measurements which are very 
difficult to make accurately. 

4. The curve of settling velocity v s for quartz grains against their 
sieve diameter (fig. 6) is taken from Rubey (17) for a temperature of 
68° F and was found to give reasonably good values. Different inves- 
tigations differ widely in the values for the settling velocity depending 
on the shape of the grains. If sediment of specific gravities different 
from that of quartz or at different temperatures is to be studied, the 
reader should refer to Rubey's original paper. 

5. The curve of £ against DfX (fig. 7) has been derived entirely 
from flume experiments with mixtures. These experiments covered 
six different mixtures but obviously did not cover all possible combi- 
nations of grain sizes. Especially, no tests were made of unsorted 

mixtures such as those found in mountain rivers near the upper end 

of alluvial stream systems where slopes are steep. This curve is 
roughly constant at £=1 for all grains with D^>X, whereas the curve 
for D<iX has a slope of about 2. This means that the lift force 
decreases about with D 2 . The same result is obtained if the effective 
velocity is assumed to decrease linearly with the size. This effect 
could be expected to occur in the laminar flow of a sublayer. Why 
the same curve seems to describe the reduction of the force in turbu- 
lent flow (X=.77A) is not clear. It may be anticipated that the 
curve for D<C,X appears to give a minimum value of transport for 
the bed under consideration; no maximum curve is known today. 
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6. The curve of Y against kj8 (fig. 8) seems to be well defined by 
flume experiments. The scatter of the points is much smaller than 
that of the £ curve. It shows the effect on the lift force as expressed 
by equation (46) if the bed as such is not hydraulicallv fully rough. 
The relationship, although entirely empirical, may be interpreted as 
a correction of the lift coefficient. 

7. The — curves (figs. 9 and 10, in pocket, inside back cover) 
are entirely theoretical, and represent equation (57). The three 
constants A*, B*, and rj 0 are obtained, yet fully supported by reliable 
experiments. For practical purposes the use of the curve in figure 
10 instead of equation (57) is satisfactory. The curve practically 
levels off at ¥*=25 even if, theoretically, <f>* becomes zero only at 

= oo . For practical purposes may be assumed to become zero 
at ^*=25. The curves in figures 9 and 10 are identical. Figure 9 
shows the comparison of the theoretical curve with flume measure- 
ments, while figure 10 is added as a work sheet with larger scales for 
easier reading. 

FLUME TESTS WITH SEDIMENT MIXTURES 

The formulas for the determination of the bed-load function, 
which cover both the bed material in suspension and that moving as 
surface creep, may be tested against flume experiments. Bed mate- 
rial goes into suspension only when it moves at very high rates and 
when its settling velocity is moderate. The experiments for checking 
the general formulas thus must be made with high-intensity flows 
over bed material of fine sand. Reports on experiments of this type, 

which incorporate all the necessary measurements, could not be 
found in the literature. As a result, the writer ran a special set of 
26 such experiments during the years 1944^46 at the Cooperative 
Laboratory of the Soil Conservation Service and the California 
Institute of Technology in Pasadena, Calif. Six different sand 
mixtures were used with various flows. 

Bed-load experiments usually are performed in a flume equipped 
with a sediment feeder at the upper end and a settling box at the 
downstream end. Sediment is added to the flow at the upstream 
end according to a predetermined rate until it is deposited at an equal 
rate downstream. In this way, an equilibrium rate is established 
for the prevailing flow conditions. In planning test experiments which 
were supposed to cover high sediment-discharge rates, such as 5 percent 
of the flow or up to about 15,000 pounds per hour, it became imme- 
diately apparent that the task of weighing, feeding, and separating 
such quantities would be far out of proportion to the scale of operations 
in a flume 10 inches wide and 30 feet long. It was concluded that 
such experiments could not be conducted feasibly unless an entirely 
different method of operation were employed. This new method 
was actually suggested by the flume itself. 

The flume was one of the so-called circulating tvpes; i. e. the water 
was driven by a propeller pump after leaving the downstream end 
of the flume. From there the water was passed through a return 
pipe directly back to the upstream end of the flume without losing 
its velocity. It was easy to operate the installation in such a way 
that the flow velocity in the entire return system could be kept higher 
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than in the flume. This gave the return system a higher capacity 
to move sediment. Actually, this capacity was so high that the 
entire load moving in the flume could be transported through the 
return channel as wash load in suspension. By this method the 
same pump actually recirculated continuously both water and 
sediment. The sediment load was measured by sampling in a vertical 
branch of the return pipe. A special study showed that the entire 
load in the return pipe was in suspension and very well distributed 
over the cross section. 

In this system, the following variables were measured: 

(1) During the run, the discharge was measured in a contraction 
of the return pipe, the slope of the water surface was measured with 
stage recorders in special pressure wells at four locations and the load 
was determined by sampling in the return pipe. 

(2) After the run, the location and slope of the bed were measured 
with a point gauge and the bed composition was obtained by sampling 
and analysis. This gave directly q T i T , S B , ib, S v and Q. Thus the 
water depth d, the average velocity v, and the slope of the energy 
grade line S e could be calculated. With the roughness of the side 
walls known, R R was calculated, and from u and Z? w of the bed the 
values of u+ and W were determined by trial and error. From these 
values and for the individual grain-size ranges were computed 
separately with the various formulas and graphs given in the pre- 
ceding chapters. 

The results, which may give the reader a conception of the reliability 
of the method, are shown in the graph of figure 11. The curve 

of that graph is the same as those in figures 9 and 10, which are given 
for comparison. A large number of points are concentrated near the 
curve whereas others scatter rather widely. In judging this scatter 
one must remember that each point represents only one sieve-size of 
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Figure 11. — graph for individual grain sizes of flume experiments. 
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one experiment. Many such grain sizes are scarce; either i B , i b or i T 
may be only a few percent. As a considerable sampling error cannot 
be avoided in the determination of the i-values, all points were omitted 
where one of the i-values was 1 percent or less. 

As the curve of equation 57, as given in figures 9, 10 and 11, 

is not used for the calculation of the individual and values, this 
graph provides additional proof that the approach which leads to 
equation 57 is basically correct. 

SAMPLE CALCULATION OF A RIVER REACH 

The greatest difficulty in applying the different equations and 
graphs to a natural river channel is the basically irregular flow in such 
a channel. Each cross section is different from all other cross sections. 
Every vertical in a cross section is different from every other vertical. 
To the writer's knowledge no usable theory exists that will permit the 

{jrediction of the flow patterns in the individual verticals. But these 
ocal flow distributions determine the local sediment motion. The 
next best approach would be a statistical description of the different 
local flow patterns, but not even that is possible with our present state 
of knowledge. We may be able to determine the statistical distribu- 
tion of the water depths, at least for existing river reaches, but the 
statistical distribution of the shear, the it^-values, which enter both 
the bed-load and the suspended-load equations, are still unknown. 

Several authors (12, 14, 18) have proposed that u+ should be con- 
sidered to be proportional to the depth d. This is equivalent to 
assuming that the local slope of the energy gradient is the same for all 
points of a channel. Even though the energy slope may have a 
tendency under certain conditions to be more constant than, for 
instance, the depth or the velocity, this assumption leads to results 
w T hich are not at all substantiated by observation. If S e were assumed 
to be constant, it would be immediately apparent from equation (6) 
that u'f must be proportional to the local W . Since the transport of 
bed load is increasing with increasing u+, the deepest points of the 
channel should always have a maximum transport. But observations 
of both flume flows and river models show distinctly that the transport 
in the points of deepest scour is usually very small or even zero and 
that the most intense transport occurs in the medium depths and on 
the top of sediment bars. This may seem to contradict a previous 
statement that the energy slope has been found to be very constant 
along bed-load-moving flume-flows. Actually, there is no con- 
tradiction. 

Although the distribution of the transport can be described at the 
different parts of a sand bar, it seems to be impossible to measure the 
energy slope locally at the different parts of any individual bar. A 
significant change of the energy level is found only between sections 
which are at least one or several bars apart. The two observations 
thus pertain to variations of a different order of magnitude. 

River observations, on the other hand, reveal that in alluvial chan- 
nels with meandering thalweg, the bars with shallow water depth 
often show a much larger transport than the deepest parts of the chan- 
nel. This can be explained only by an irregular distribution of the 
energy slope, both in the cross section and along the flow lines. For 
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want of better information about the distribution of the shear stress 
over, the entire bed area, therefore, it is proposed to use the average 
flow conditions in the description of the channel behavior. Except 
in extremely wide and flat channels, such as on debris cones, this 
method seems to give reasonably good results. 

Choice of a River Reach 

In practical calculations of the bed-load function for a particular 
river reach, the length of the reach must be sufficient to permit ade- 
quate definition of the over-all slope of the channel. The channel 
itself should be sufficiently uniform in shape, sediment composition, 
slope and outside effects such as vegetation on the banks and over- 
banks, that it can be treated as a uniform channel characterized by an 
over-all slope and by an average representative cross section. Such 
a section can be described by two curves in which the cross-sectional 
area and the hydraulic radius are plotted against the stage. 

Description of a River Reach 

One problem is that of determining how a number of cross sections 
can best be averaged. As the river reach is to be treated as a uniform 
channel with Constant cross section and slope, in which only uniform 
flows are studied, a representative or average slope must be found, 
together with the average section. If a sufficiently long and regular 
profile exists for the river under consideration, the general slope of the 
reach should be taken from it. In the absence of such a profile the 
slope must be derived from the cross sections themselves. Under all 
conditions, the cross sections must be tied together by a traverse 
which gives their relative elevations and the distance between them 
along the stream axis. Then the wetted perimeter and the wetted 
area are calculated for various water surface elevations. These are 
plotted in terms of the water surface elevation for each cross section. 

It is fairly common usage to construct the stream profile from the 
lowest points of the sections. This procedure is satisfactory for a 
long profile. If the reach is short, however, the use of a low-water 
surface is more satisfactory as the influence of insignificant local scour- 
holes is excluded. If such a low-water profile is not recorded when the 
sections are surveyed, a profile found from the area-curves may be 
substituted. A characteristic low-water discharge may be selected 
for the streams. The average velocity for such a flow can be estimated 
roughly. By division of the two one may find the corresponding 
low-water area of the cross sections. If the water-surface points 
which give this area at the different sections are connected, an approxi- 
mate low-water surface is defined which represents a profile that is 
more regular and more representative than the profile of the low points 
of the bed. 

After the representative slope is selected by fitting a straight line 
through the profile points, this slope may be used in averaging the 
cross sections. This can be done by sliding all the sections along this 
average slope line together into, for instance, the lowest section. 
With the sections described as two curves, namely (1) of the area, and 
(2) of the wetted perimeter, both in terms of water-surface elevation 
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and all reduced into one plane, the areas and the wetted perimeters 
for each elevation may then be averaged directly. By this means, 
averages are obtained for areas and wetted perimeters which corre- 
spond to a water-surface line inclined according to the representative 
slope. This makes the procedure consistent with the assumption of 
uniform flow. 

The averaged area and perimeter curves allow the direct calculation 
of the hydraulic radius R in terms of the stage. The average flow 
velocity is then calculated for these R values and the discharge is 
obtained by multiplication with the area. Each R value is thus 
assumed to be representative for a discharge over the entire reach. 

Application of Procedure to Big Sand Creek, Miss. 

The procedure as outlined so far has been applied to a reach of 
Big Sand Creek, a notorious sediment carrier near Greenwood, Miss. 
This stream has a characteristic fine-sand alluvial bed although it 
drains only about 100 square miles. 

Hydraulic calculations for channels in general 

Step 1: The location of the reach was determined by the location 
of the only existing gaging station on Big Sand Creek. The length 
selected, about 3 miles, was based upon the following considerations: 
The reach must be treated as a uniform channel with uniform flow; 
it must be possible, therefore, to neglect any changes of the velocity 
head u 2 /2g affected by the total energy drop in the reach. At high 
water, the flow velocity is 10 to 11 feet per second, for which value 
the velocity head is about 2 feet. A 30-percent variation of the 
velocity head, or about 0.6 feet, may be expected to occur within the 
reach. The total energy drop in the reach with a length of 3 miles 
and a slope of 0.001 is 15 feet. The uncertainty of the effective slope 
is thus 0.6/15=0.04 or 4 percent, deemed to be tolerable. 

Step 2: Cross sections, well distributed over the 3-mile reach, were 
surveyed and plotted (fig. 12). The distances between sections are 
shown in the profile (fig. 13). The elevations plotted at all sections 
of figure 13 refer to the low-water level with a wetted area of 50 square 
feet. The value 50 was chosen as it is representative of a low-water 
flow in this stream. The elevations themselves are taken from figure 
14, which gives for all cross sections the wetted area in terms of water- 
surface elevation. These values may be derived from the cross sec- 
tions of figure 12. Similarly, from the same cross sections, a graph 
may be derived which gives the wetted perimeter in terms of water- 
surface elevation. In the case of this wide and shallow channel, the 
wetted perimeter was assumed to equal the surface width. 

Step 8: The most probable straight line was then laid through the 
50-square-foot points of the profile shown in figure 13 by the method 
of least squares. The slope of this line was determined to be 5.54 
feet per mile or 0.00105. This represents the slope of the uniform 
channel which will be used in the determination of the bed-load 
function. 

Step 4* The average or representative cross section was determined 
by sliding all cross sections down the channel along the slope 
£=0.00105 into the plane of the section at the lower end of the reach 
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Figure 12.— Cross sections of Big Sand Creek, Miss., in their actual position. 
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Distance along C of the channel (1,000 feet) 

Figure 13— Profile of Big Sand Creek, Miss., showing elevations of water 
surface at cross sections 40 to 51 for which the wetted area is 50 square feet. 
Straight-line profile obtained by method of least squares. 
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Figure 14. — Area stage curves for sections 40 to 51 of Big Sand Creek, Miss., 

in their actual position. 

(fig. 15). This was accomplished by lowering each curve of figure 14 
by the distance 0.00105 L x in which L z is the distance of the cross 
section from the downstream end of the reach measured along the 
stream axis. Figure 15 gives these transposed sections and the aver- 
age section for all stages. The heavy curve connects the averages 
and represents the area curve of the average or representative cross 
section. In a similar way the curve of the average wetted perimeter 
in terms of stage was determined and combined with the area curve 
(fig. 16). Division of corresponding values from the two curves 
makes it possible to determine, point by point, the curve of the 
hydraulic radius against the stage (fig. 16). 

It is usually impossible to construct a cross section which satisfies 
both the curves of area against stage and of hydraulic radius against 
stage. The procedure illustrated in figure 16 is fully satisfactory, 
however, as it develops the most representative average section with 
respect to its hydraulic behavior. Where no outside roughness 
such as islands, submerged rocks, or rough banks exist in the channel, 
these two curves are sufficient to define the shape and slope character- 
istics of the channel. 

Step 5: The grain-size composition of the bed is determined by 
sampling. A bed which appears to be very uniform, such as that of 
Big Sand Creek, may be described by three to five samples. Each 
of the four samples listed in table 5 was a composite of three or four 
cores, taken in tne same cross section at evenly spaced points over the 
total width of the channel. The individual samples were obtained 
by means of an auger or a pipe-sampler and were taken down to 
a depth of about 2 feet, the estimated depth of scour or active bed 
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Figure 16. — Description of the average cross section, Big Sand Creek, Miss. 
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Figure 17— Average-grain-size analysis of the bed, Big Sand Creek, Miss. 
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movement. The four samples referred to in table 5 thus represent 
the average grain-size composition at four cross sections which them- 
selves are well distributed over the length of the reach under consider- 
ation. From figure 17, which is a logarithmic-probability plot of 
the analyses of the data shown in table 5, the characteristic grain 
sizes of the bed may be read. The size which enters the equations 
of transport is D 35 =0.29 millimeters =0.00094 feet (35 percent of the 
mixture is finer); and the size characteristic for friction D 65 =0.35 
millimeters =0.001 15 feet (65 percent of the mixture is finer). 

For the description of more heterogeneous beds, more and larger- 
sized samples are required. For gravel-carrying beds it may be 
necessary to analyze samples of several hundred pounds, especially 
if the spread of grain sizes is large. 



Table 5. — Distribution of grain sizes in 4 bed samples, Big Sand 

Creek, Miss. 



Grain size in 
millimeters 


Sample 
1 


Sample 
2 


Sample 
3 


Sample 
4 


Average 














Milli- 




Percent 


Percent 


Percent 


Percent 


Percent 


meters 


D> 0.589 


1. 3 


1. 4 


5.7 


1. 4 


2. 4 




0.589>D> .417 


18. 6 


9. 8 


27. 5 


15. 4 


17. 8 


0. 495 


.417>D> .295 


47. 7 


36. 1 


35. 9 


40. 9 


40. 2 


. 351 


.295>D> .208 


28. 2 


40. 4 


23. 7 


35. 6 


32. 0 


. 248 


.208>D> .147 


2. 6 


9. 3 


5. 7 


5. 7 


5. 8 


. 175 


.147>D 


1. 6 


3. 0 


1. 5 


1. 0 


1. 8 





Step 6: As the sediment transport will be calculated for the grain 
sizes between 0.589 millimeters and 0.147 millimeters, which cover 
95.8 percent of the bed material, the calculation will be made for indi- 
vidual sieve fractions using as representative the average grain sizes 
of 0.495, 0.351, 0.248 and 0.175 millimeters or 0.00162, 0.00115, 
0.00080 and 0.00057 feet, respectively. 

Hydraulic calculations for channel without bank friction 

Step 7: An analysis of the equations shows that the most direct 
approach is obtained if values of the hydraulic radius with respect 
to the grain, R b ' , are assumed (column 1, table 6). 

Step 8: From R b ' the corresponding friction velocity u* is cal- 
cidated using equation (6), page 9. 4 

Step 9: The thickness of the laminar sublayer 5 is obtained, in feet, 
from equation (5), page 8. The kinematic viscosity v is assumed to 
have the value 10~ b ft. 2 /sec, which is correct at a temperature of 75° 
F. The friction velocity u+ is taken from column 2, table 6. 

Step 10: With £55= 0.00 115 feet, the values of kJ8 (column 4, 
table 6) are calculated. 

4 Equation (6) as well as all other equations in this publication are dimen- 
sionally correct so that any consistent set of units may be used. In table 6 and 
all following tables, the engineering system of foot-second-pound units is used. 
The value of ii* thus is obtained in feet per second (ft./sec), if S t is introduced as 
a tangent (£.=0.00105), R b -' in feet, and 0=32.2 ft./sec*. 
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Step 11: The correction x for the transition from smooth to rough 
boundaries may be read from the graph of figure 4 (column 5, table 6). 

Step 12: The values of the apparent roughness A=k,/x are calcu- 
lated in feet (column 6, table 6). 

Step IS: Now, the average flow velocity u is calculated from 
equation (9), page 10 in feet per second (column 7, table 6). 

Step 14' Next, for the determination of the frictional contribution 
of the channel irregularities, the parameter S^' is calculated according 
to equation (11), page 10 (column 8, table 6). The parameter is 
dimensionless, as is the ratio of the densities, which may be determined 
as follows for Big Sand Creek: 

s l=Il = 'JlLZl = s, - 1 - 2 . 65 - 1 - 1 . 65 
s / 1 

in which S, is the specific gravity of the solids and has a value between 
2.65 and 2.68 for most natural sediments. D n and R b must be 
entered in the same units, for instance in feet (Z> 35 = 0.00094 feet), and 
the slope S e of the energy gradient as an absolute figure (5^=0.00105 
for Big Sand Creek) . 

Step 15: From figure 2, the values of u/u+" (column 9, table 6) are 
read for these values. 

Step 16: uj' is calculated in feet per second (column 10, table 6). 

Step 17: Equation (6), page 9, in the form 

allows the calculation of R* (column 11, table 6) in feet from u" 
(column 10, table 6), £,=0.00105 and 0=32.2 ft./sec. 2 . 

Step 18: The two components R b and R" are usually the only 
components of the hydraulic radius R b pertaining to the bed and may 
be added directly 

R b =R b '+R b " (64) 

as is done in column 12, table 6. All R values are measured in feet. 

Step 19: Where no additional friction such as that from banks or 
vegetation must be introduced, R b represents the total R of the 
section. In this case, the total cross-sectional area A T in square feet 
(column 14, table 6) and the channel width (wetted perimeter) p b in 
feet (column 15, table 6) are read directly from figure 16 for each 
R value. 

Step 20: The flow discharge Q in cubic feet per second (column 16, 
table 6) is calculated as 

Q = A T -u (65) 

Step 21: The rating curve of figure 18 for the average section, as 
shown in figure 16, is obtained by plotting the discharge Q (column 16, 
table 6) against the stage (column 13, table 6), which itself is read 
directly from figure 16 as a function of R. As the entire discharge 
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rating curve is needed, it is not too important which points are 
actually calculated. The choice of the points may be made, therefore, 
according to the greatest ease of calculation. 




100,000 



Figure 18. — Rating curve of the average cross section, Big Sand Creek, Miss. 

Hydraulic calculations for channel with bank friction 

Whenever a channel has wetted boundary areas which consist of 
material different from the movable bed material, or if some wetted 
boundary areas are covered with permanent vegetation, these areas 
represent what have been called "bank surfaces" and must be intro- 
duced separately in the calculation. Whether a given small percent- 
age of "bank surface" has an appreciable influence, or not, can be 
determined only by a trial calculation. 

The calculation with bank friction is usually somewhat complicated 
by the fact that this additional bank friction must be considered in 
terms of the stage and not in terms of the bed friction R b ' or R b . A 
trial-and-error method must therefore be used for its solution. Up 
to and including the determination of R b by equation (64) (step 18), 
this calculation is identical to that without bank friction; but with 
bank friction, R b is not equal to the total hvdraulic radius R of the 
section. Instead of this simple equality, the procedure previously 
outlined by the author (8) must be used. The banks are assigned a 
separate part A w of the total cross section A T , a wetted perimeter p M 
of the bank surface and a hydraulic radius R w defined by equation (66). 

A m =R w • p w (66) 
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If A b is that part of the cross-sectional area pertaining to the bed and 
if no friction acts on the flow except that on the bed and that on the 
banks, the following equation holds: 

A T =A b +A„ (67) 

in which A T is the total area of the cross section. The partial areas 
may be expressed by the hydraulic radii 

A T =p b Rb+pw R* (68) 

The hydraulic radius of the bed R b is calculated in terms of R b ' (column 
12, table 6) while R w may be calculated for each R b ' value and the 
corresponding average velocity u from equation (69) 

for any chosen roughness of the banks n w . This assumes that the 
average flow velocity in all parts of the cross section is the same as 
that of the total section and that the proper friction formula may be 
applied to each part of the cross section according to its frictional 
surface independent of the friction conditions in the remainder of the 
cross section. Choosing, for instance, a value of n tp =0.050 for the 
banks and using the average energy slope <S=0.00105 and the average 
velocities of column 7, table 6, the R„ values of column 17 result, also 
in function of R b ' . 

The p values of equation (68) are not determined in function of R b ' f 
however, but in terms of the stage of figure 16. The curve of p b has 
been adequately described. The curve for p w is assumed arbitrarily 
as a straight line defining p w as twice the height of the water surface 
above elevation 150.0 feet. The Big Sand Creek cross sections do not 
give information on the bank roughness, but photographs and some 
notes made by a field inspection gave a basis for developing this curve, 
as well as the assumed value for n w =0.050. 

Alternate step 18: With p b and p w functions of the stage, which 
itself is related in figure 16 to the total area A T , with R b and R w a func- 
tion of Rb as the flow velocity, and with equation (68) tying them to 
A T and the p values, only a trial-and-error method may be used to 
find solutions. Figure 19 shows how this is done. The A T stage-curve 
of figure 16 is intersected by short curves, the points of which are 
obtained by calculating A T values according to equation (68), assuming 
for each such curve the R b and R w values according to one of the cal- 
culated Rb values, with p b and p w varying according to the stage. 
Plotting the calculated (p b Rb+Vw Rw) values against the assumed 
stages, it is clear that the actual A T points are found for the calculated 
Rb values at the intersection points with the A T curve. 

Alternate step 19: The elevations of intersection are given in figure 
19 and used in column 18 of table 6 to determine Pb (column 19, table 
6) and A T (column 20, table 6). 

Alternate step 20: The discharge Q (column 21, table 6) is calculated, 
using equation (65), page 53. 
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Figure 19.— Trial-and-error determination of the rating curve with bank friction 

Big Sand Creek, Miss. 

Alternate step 21 : The result is plotted in figure 18 together with the 
curve for the same section neglecting bank friction. It is interesting 
to note that the discharge for a given R b changes very distinctly with 
the introduction of bank friction, but that the rating curve is not 
greatly affected. The values of R Tt which must be used as the average 
depth in the calculation of A for the suspension-integrals, are defined 
arbitrarily by 

R T =A T lp b (70) 

and are given in column 22, table 6. This completes the purely 
hydraulic calculations of the river reach with bank friction. 

Sediment-rate calculation 

The sediment transport is calculated for the individual grain-size 
fractions of the bed and for the entire range of discharges. In this 
connection, it is advantageous to distinguish in a separate table the 
steps which are common to all grain sizes from those which must be 
performed separately for each grain size. As the calculation of trans- 

{>ort rates is performed for the flow rates given in table 6, these calcu- 
ations are added to that table as columns 23 to 27. The transport 
calculations for the channel without bank friction are shown in table 
7, and those with bank friction appear in table 8. 

These two calculations are not significantly different; e. g., for the 
determination of A without bank friction the hydraulic radius R b is 
used (table 7), whereas R T is used when bank friction is a factor 
(table 8). As comparison of the two tables will show, the only 
important differences between the two calculations are in the strictly 
hydraulic relationships, as in the variation of p b and of Q for equal 
R b values. 
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The sequence of operations and the application of formulas for the 
calculation of tables 7 and 8 are as follows: 

Step 22: The representative grain sizes are chosen at the geometric 
mean of the two size limits of each sieve fraction (column 1, tables 
7 and 8). 

Step 23: The fraction i b of the bed material for the same sieve 
fraction is determined as an average of all bed samples available for 
the river reach in question (column 2, tables 7 and 8). 

Step 24: The same R b values are used as in table 6 (column 3, 
tables 7 and 8). 

Step 25: Next, ¥ is calculated using equation (49) page 36 and 
Sj -^=S s —l, as demonstrated previously (column 4, tables 7 and 8). 

Step 26: The characteristic distance X is derived according to the 
relationships (45), page 35, in feet, with A and 8 given in columns 6 
and 3 of table 6, respectively (column 23, table 6). 

Step 27: The ratio D/X is calculated from D in column 1 and X in 
column 23, both from table 6 (column 5, tables 7 and 8). 

Step 28: The £ values are read from figure 7 in terms of D/X (column 
6, tables 7 and 8). 

Step 29: The Y values are read from figure 8 in terms of which 
is found in column 4 of table 6 (column 24, table 6). 

Step 30: fi x is calculated from equation (49), page 36, with A" given 
in column 23 and A in column 6, both from table 6 (column 25, table 6). 

Step 31: j8=log 10 (10.6) = 1.025 is divided by /3 X of column 25 and 
the result is squared (column 26, table 6). 

Step 32: The parameter ^ is calculated according to equation 
(54), page 37, with ^ in column 4, table 6; £ in column 6, tables 7 and 8; 
Y in column 24, tables 7 and 8, and (/3/ftj-) 2 in column 26, table 6 
(column 7, tables 7 and 8). 

Step 33: From figure 10, is read in function of SF* (column 8, 
tables 7 and 8). 

Step 34: The bed-load rate, i B q B , is calculated from using 
equation (41) and (42), page 34, in the form 

where: 

<P+ is given in column 8, tables 7 and 8 ; 
i b in column 2, tables 7 and 8 ; 

s, is from the mechanical analysis, but may usually be assumed 

to be 2.67-1.93 = 5.17 slugs/ft. 3 ; 
0=32.2 ft./sec. 2 ; 

D is from column 1, tables 7 and 8; and 
(S s — 1)* is 1.29 (column 9, tables 7 and 8). 

Step 35: A is calculated as A=2D/R b in table 7 and as A=2D/R T 
in table 8 (column 10, tables 7 and 8). 

Step 36: The exponent z is calculated according to equation (27), 
page 17, where v t may be read from figure 6 for various grain sizes. 
These curves apply to average-shaped quartz grains at a temperature 
of 20°C. For other materials and temperatures the settling velocities 
may be measured or calculated from Rubey's formula (17). The 
constant 0.40 seems to be a usable average (19), while u*' is given in 
column 2, table 6 (column 11, tables 7 and 8). 
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Step 87: The integral I x is read from the diagram in figure 1 for 
the various values of A and z (column 12, tables 7 and 8). 

Step 88: The integral I 2 is read from the diagram in figure 2 for 
various values of A and z (column 13, tables 7 and 8). 

Step 89: The coefficient P of I x is calculated according to equation 
(62), page 40, where: 

x is from column 5, table 6 ; 
k t =Dn= 0.001 15 feet; and 

d equals K b in table 7 and R T in table 8 (column 27, table 6). 

Step 40: The expression (P/1+/2 + I) is calculated according to 
equation 63, page 40, (column 14, tables 7 and 8). 

Step 41 : The total transport rate per unit width and time for the 
individual size fraction i T q T is calculated according to equation 63, 
page 40, in lbs. /ft. -sec, (column 15, tables 7 and 8). 

Step 4%: The same total transport rate reduced to r&=l, which is 
helpful in judging the behavior of the different grain sizes of a bed 
mixture ir^r/ib, is calculated in column 16, tables 7 and 8. 

Step 48: The total sediment transport rate for the entire section 
is calculated for individual size fractions in tons per day from the 
equation 

i T Q r =(irq T ) 43.1 p b 

where i^r appears in column 15, tables 7 and 8; and p b in column 15 
of table 6 if bank friction is neglected, and in column 19, table 6 if 
bank friction is introduced (column 17, tables 7 and 8). 

Step 44'- Finally, the total transport rates are added for all sedi- 
ment sizes coarser than any size Z), in tons per day (column 18, tables 
7 and 8). This represents the final form in which the bed-load func- 
tion is presented (see figs. 22 and 23). 

Discussion of Calculations 

The behavior of the individual grain sizes within the sediment 
mixture is best characterized by the curves of figure 20 for Big Sand 

Creek. They give the values %T '. gr against grain size for various 

stages. Each curve thus refers to one discharge and gives the rate 
at which the individual grain sizes would move if they covered indi- 
vidually the entire bed area. The size range has been extended far 
beyond the range that is actually important in Big Sand Creek in 
order to show the characteristic parts of these curves. 

Figure 21 gives comparative curves for the Missouri River. The 
graphs for both streams are calculated without bank friction. It is 
apparent that both sets of curves have a very similar character 
although the curves themselves are distinctly different; viz, they 
both show one well-defined maximum and a tendency toward a second. 
Because of their rather characteristic shape such curves have been 
called camel-back curves (c. b. c). 

It appears that only few grain sizes, D, must be calculated to define 
these entire curves, from which intermediate points may then be 
interpolated. 
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Figure 20. — Curves of against D for various stages, Big Sand Creek, Miss. 

Regardless of whether or not the entire range of sizes has practical 
importance, it is useful to evaluate the significance of the different 
parts of the curves. At the maximum values (about 0.1 millimeter) 
the material moves almost exclusively in suspension. The values of 
i T q-r/iB Q_b are far above 1,000 for the higher stages for which this 
maximum is especially pronounced. The drop to the left of this maxi- 
mum becomes very steep and is caused by the relatively fast increase 
of £ and with it of with decreasing D. This sudden drop of the 
transport occurs after increases above 20 as shown in figure 10. 
A similarly steep drop occurs at the right end of the curves: above 
D=10 millimeters for Big Sand Creek and around 100 millimeters 
for the Missouri River at flood stages. 

The reason for this drop is the increase of above 20 due to growing 
D values. The transport curve thus is limited at both ends of the 
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Figure 21. — Curves of against D for various Missouri River stages. 
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D range by increased values. But the reason for this increase 
at the two ends is different. At large diameters, the weight of the 
particles becomes too great for the lift forces, whereas at small diame- 
ters the effective velocities of flow and turbulence are reduced as the 
grains begin to hide— that is, come to rest beyond the influence of 
turbulence between larger grains or in the laminar sublayer. 

The usually very steep increase of irQrlib from D=l millimeter 
down to 0.1 millimeter, especially at higher stages, is caused by an 
important shift from surface creep to suspension. This transition 
occurs with a change of the exponent z from about 2 to 0.2. There 
remains only to explain the decrease of ir^r/ib between a size of about 
10 millimeters (see Missouri River curve, figure 21) down to about 
1 millimeter. The tendency to decrease in this range is again espe- 
cially pronounced at flood stages. This reduction is still in the range 
where suspension is unimportant. It must be explained by bed-load 
motion only. It results from the condition where the frequency * 
of individual motions of particles increases more slowly with decreas- 
ing D than the volume of the individual particles and the length of 
their jump is reduced. In other words, larger numbers of the smaller 
particles move but the rate of their movement by weight remains 
smaller because the individual particle weight decreases very rapidly 
with D. Commonly, for instance, for all curves of Big Sand Creek 
and for the low stages of the Missouri River, this intermediate mini- 
mum or dip in the camel-back curve is not pronounced because the 
beginning of increasing suspension overlaps with the decrease of bed- 
load rates due to decreasing particle size. Thus, the detail shape of 
the camel-back curves varies considerably as well as the location of 
maxima. 

How may the rather large rates of transport in the silt and clay 
sizes of many streams be explained in the light of this analysis? The 
camel-back curves for many streams decline just as sharply near 0.01 
millimeter as do the two examples. This becomes more understand- 
able as one recalls that the method of calculating sediment loads 
presented in this publication tends to give minimum rates, as did 
the experiments from which the method is derived. This is reflected 
in the £ curve, which illustrates the fact that the small bed-particles 
always do their utmost to-hide behind larger grains or in the laminar- 
sublayer. Thus, it is easy to visualize that a rather small amount of 
similar particles in more prominent positions of the bed could support 
a very large additional transport of these particles without increasing 
effectively their over-all concentration i b in the bed as a whole. 
This might explain the existence of wash load on the same bed which 
sustains a bed-load function. No quantitative information is yet 
available on this condition, however. 

From the camel-back curves, or directly from the sediment load 
calculations, the actual transportation rates ir Q r of the entire cross 
sections may be determined for the calculated flows. These values 
i T Q t are calculated for individual size fractions which themselves are 
chosen arbitrarily. As in the case of mechanical analyses, the use of 
cumulative curves is more general than that of distribution curves. 
The information shown in columns 17 and 18, therefore, is given in 
figures 22 and 23 as the sum of the transport of all grain sizes coarser 
than a given limit in terms of the discharge Q. Each curve corre- 
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sponds, therefore, to an individual point of a cumulative curve. The 
group of curves gives, then, the full description of the bed-load 
function. 

The fact that the curves, especially the curves for the total trans- 
port of all particles coarser than 0.147 millimeters, are much flatter 
than 45° indicates that the sediment concentration increases with 
rising stage. Calculating the load in parts per million (p. p. m.), the 
channel without bank friction gives, for instance, at 400 cubic feet 
per second, a concentration c of 

670 • 10 6 - 2 000 Ano 
C= 49C 62.4- 3600- 24 = 6 ° 9 Ppm ' 

The corresponding concentration at a flood stage of 30,000 cubic feet 
per second is 23,900 parts per million. Each cubic foot of water at 
highest flood stage is, therefore, about 40 times as effective in moving 
sediment as it is at a rather low stage of 400 cubic feet per second. 
In order to predict the load that will be moved by a stream, it is thus 
necessary to know not only how much flow will occur but the duration 
of each rate of flow. Fortunately, the sequence in which the various 
flows occur is not important as long as an equilibrium condition is 
assumed for all flows. The most advantageous description of the 
flow-conditions is in this case the well known flow-duration curve. 

More work needs to be done in the development and publication of 
flow-duration curves for both larger and smaller rivers. Little is 
known, in particular, about the relationship between the flow-duration 
curves for tributaries of different size in the same watershed, between 
those of different sections of the same river and sections with different- 
sized contributing watershed areas: nor about the characteristics of 
the curves in different parts of the country. There is a great need for 
further development in this almost untouched aspect of hydrology (4) . 

Because most existing flow records are given as daily averages, 
tons per day was chosen as the unit in columns 17 and i8 of tables 
7 and 8. New and perhaps better methods of integrating the bed- 
load function with flow data for application to different types of 
channels may be apparent to the practicing engineer after he has 
acquainted himself with the basic concepts presented in this publication. 

LIMITATIONS OF THE METHOD 

Not only the possibilities for application, but also the limitations 
and cautions to be observed in applying any newly developed engineer- 
ing method should be presented as objectively as possible by the 
author of that method. Computation of the bed-load function, as 
presented in this publication, undoubtedly can be and will be improved 
and its application extended as the results of additional basic research 
become available and as the method is applied to a greater range 
of practical field problems. It seems desirable here, however, to 
point out certain limitations that can now be recognized and which 
may serve as a stimulus to further research and field trials needed 
to remove these limitations. 

In an effort to devise a unified method of calculating the transport 
of bed sediment for immediate practical application, it was necessary 



Digitized by Google 



66 TECHNICAL BULLETIN 1026, U. S. DEPT. OF AGRICULTURE 



to develop several strictly empirical relationships such as the curves 
of figures 5, 6, 7, and 8, for which only a very limited amount of 
substantiating data now exist. Other, somewhat different but 
practically equivalent systems could be used in developing these 
curves and in defining X. Preliminary calculations indicate that 
such alternate systems give essentially similar results when applied 
to large river conditions. The author has chosen the particular 
system of curves and parameters used herein because they seem to 
have a greater significance within the framework of fluid dynamics 
theory. It is recognized, however, that a better set of curves may 
be, and probably will be developed as more information becomes avail- 
able. Inasmuch as these curves give results that appear to be con- 
firmed by all the field checks thus far made, they do not seem to be 
open to criticism so much from the standpoint of accuracy as from 
lack of a solid theoretical foundation, such as exists for the integrals 
of suspension, the <p* — curve, and the hydraulic curve of figure 4. 

The justification for this publication, other than the urgent need 
for a usable method of computing bed-load transport, rests, in the 
author's opinion, on the basic soundness of the two major principles. 
These are (1) the restriction of bed-load relationships to the bed layer, 
and (2) the method of relating the movement of bed material in 
suspension to the concentration in the bed layer. Even though some 
of the constants and approximations used in defining these relation- 
ships may be subject to later improvement and refinement, the 
unified method as presented appears to be basically correct. It 
has been confirmed in all cases where it has been checked against 
the load transported through a river reach and subsequently deposited 
in a place where its volume could be measured. 

It must be clearly understood that this method does not permit 
the calculation of the total sediment load of a stream, but only the 
total transport in suspension and in the bed layer of bed-material 
load through an alluvial bed channel. That part of the total load 
which is not included in the bed-load function cannot be determined 
by any analytical method now known. It must be measured by 
suspended-load sampling techniques, which simultaneously measure 
that part of the bed material load that is moving in suspension at 
the time of sampling. The two parts of the load moving in suspen- 
sion can be calculated separately, however, if size analyses of the 
suspended-load samples are made. The finer fractions of the total 
sediment load moving in suspension — the wash load — do not appear 
to be a function of the flow, except in a veiy vague form. Their rate 
of transport is related primarily to the supply available from water- 
shed lands, stream banks, etc. The stream's capacity for transporting 
the finer fractions of the total load is nearly always vastly in excess 
of the supply available to it and therefore the two cannot be func- 
tionally related from a practical standpoint. 

It might appear that the method presented herein is too complicated 
and time-consuming for practical use. It is conceded that the calcu- 
lations take more time and understanding than those required, for 
example, to compute the water discharge of the stream. On the 
other hand they are no more extensive than the calculations generally 
required for a single bridge spanning a river, and are generally less 
costly than sampling the suspended load of a river at a single cross 

I 
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section for a year. Moreover, the complete sequence of calculations, 
designed for application to all sizes of streams and to the complete 
range of sediment sizes from gravel through fine sand, may be simpli- 
fied considerably for any given stream after it becomes apparent that 
some aspects of the problem do not apply. 

Lastly, it must be recognized that a certain amount of practical 
experience with, or understanding of, river behavior is highly valuable 
in the correct application of this or any other method of sediment- 
load determination. Sediment movement and river behavior are 
inherently complex natural phenomena involving a great many 
variables. The solution of practical problems cannot be simplified 
beyond a certain point. For example, judgment based on experience 
must be applied in choosing river reaches to be calculated; in deciding 
whether a given reach is actually alluvial in character and thus has 
a well defined bed-load function; in deciding whether an apparently 
alluvial stream becomes nonalluvial in character at flood stage such 
that all of the bed material moves in suspension leaving a clean rock 
bed; in the choice of some constants; etc. 

Much of the practical experience needed in any study of river be- 
havior can come only from prolonged and careful study of rivers and 
of the data obtained from them in the field. Preferably such study 
should be under the guidance of one of the relatively few experienced 
engineers now engaged in this field of activity. It seems likely that 
as more experience is gained by more engineers in the application of 
this or other methods, the experience can be systematized into forms 
that will facilitate training and understanding. For the next few 
years, at least, river problems will continue to tax the ingenuity of 
even the most highly versed specialists in this field. 

SUMMARY 

(1) A unified method of calculating the part of the sediment load in 
an alluvial stream that is responsible for maintaining the channel in 
equilibrium, namely, the bed-material load, is set forth. 

(2) The relationship between the rate of transport of bed-material 
load, its size composition, and the flow discharge is called the bed-load 
function and is explained for the case of a channel in equilibrium. 

(3) The first part of the calculation covers the hydraulic description 
of the flow for each discharge. 

(4) The resulting equilibrium transport is divided into two parts: 
(a) the suspended load which includes all particles the weight of which 
is supported by the fluid flow, and which has beon found to include all 
particles moving two diameters above the bed or higher; and (b) the 
bed load which includes all particles moving in the bed layer, a layer 
two diameters thick along the bed. The weight of all particles moving 
in the bed layer is supported by the bed as they are rolling or sliding 
along. On the basis 01 this definition, the thickness of the bed layer 
is different for the various grain sizes of a sediment mixture. 

(5) The motion of bed-material load in suspension is described by 
the commonly accepted method based on the exchange theory of 
turbulent flow. The transport is integrated over a vertical. 

(6) The description of the bed-load motion in the bed layer is the 
same for fine sand as for coarse particles which never go into suspen- 
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sion. The effect of varying ratios between the grain size and the lami- 
nar sublayer thickness must be allowed for and evaluated, however. 

(7) A complete sample calculation for a reach of Big Sand Creek, 
Miss., demonstrates the practical application of the method and of its 
formulas and graphs. 
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APPENDIX 
List of Symbols 

The symbols used in this publication, together with the page on 

which each symbol is defined and the equation in which the symbol 
is of major importance, are given in the following list. 

Equation 

Page No. 

0 Thickness of bed layer 17 (28) 

A=afd Ratio of bed-layer thickness to water depth (dimen- 

sionless integration limit of suspension) 18 (33) 

A\ Constant of grain area 33 

A 2 Constant of grain volume 32 

A 3 Constant of time scale 33 (37) 

A s Constant of bed-layer concentration 39 

A 6 Constant of bed-layer concentration 40 

A i Cross-sectional area pertaining to bed 55 

A c Cross-sectional area, synonymous with A T .- _ 55 (67) 

A L Constant of the bed-load unit-step__ 33 (38) 

A T Total area of a cross section 55 (67) 

A w Part of the cross section pertaining to the banks 54 (66) 

A* Constant, scale of 34 (41) 

A' Cross-sectional area pertaining to the grain 9 

A" Cross-sectional area pertaining to irregularities 9 

B Constant, scale of * 36 (49) 

B' Constant, scale of * 36 (51) 

£* Constant, scale of** 37 (54) 

c Concentration in dry weight per unit of volume 15 

c 0 Concentration at distance a from bed 17 (29) 

cl Lift coefficient 31 

c„ Concentration at distance y from bed 17 

d Water depth 16 

D Grain size; diameter of balls 5 

D35 Grain size of which 35 percent is finer 10 (11) 

D K Grain size of which 65 percent is finer 8 

g Acceleration due to gravity 8 

i b Fraction of bed material in a given grain size 33 

i B Fraction of bed load in a given grain size- 32 

i. Fraction of suspension in a given grain size 40 (61) 

t'r Fraction of total load in a given grain size 40 (63) 

h Integral value 24 (35) 

I, Integral value 24 (35) 

A;, Roughness diameter 8 (2) 

1 A distance 4 

L Lift force on bed particle 35 (44) 

l e Distance of exchange, mixing length 15 

L z A distance in direction of the flow 48 

n w Friction factor (Manning) of the banks 55 (69) 

p Probability of a grain to be eroded 33 

P Parameter of total transport 40 (62) 

p b Wetted perimeter of the bed 9 

p L A lift pressure 31 (36) 

p. Probability of a grain to be eroded per second 33 

p w Wetted perimeter of the banks 54 (66) 

Q Flow discharge :_. 53 (65) 

qa Bed-load rate in weight per unit of time and width 32 

q, Corresponding suspended load rate 18 (31) 

q T Corresponding total load rate 40 (63) 

Qr Total sediment load in cross section 60 

3. Vertical exchange discharge per unit area 15 
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Equation 
Page No. 

R , Hydraulic radius 8 

R' Hydraulic radius with respect to the grain 9 (6) 

R" Hydraulic radius for channel irregularities 9 (6) 

R b Hydraulic radius with respect to the bed 53 (64) 

R T Hydraulic radius of the total section 56 (70) 

J?„ Hydraulic radius with respect to the bank 54 (66) 

8 Slope 46 

S B Bed slope 43 

S, Energy slope 8 

S 9 Slope of the water surface 43 

S, The specific gravity of solids 53 

t Variable of integration _ 37 (56) 

It Exchange time of bed-load particles 33 (37) 

u Velocity in direction of the main flow. _ 9 

u Time average of u, averaged over the vertical 10 

Ha Velocity with which bed load moves 39 

Hy Time average of the velocity u at y above the bed 17 (30) 

u. Shear velocity 8 (4) 

u'. Shear velocity with respect to the grain 9 (6) 

u"» Shear velocity for channel irregularities 9 (6) 

M» Velocity at the edge of the laminar sublayer _ 11 (14) 

v A vertical velocity __ _ 15 

v, Settling velocity of a sediment particle 15 

W Weight of sediment particle under water 35 (43) 

x Parameter for transition smooth-rough 8 (3) 

X Characteristic grain size of mixture 35 (45) 

y Distance above the bed 8 

Y Pressure correction in transition smooth-rough 36 (51) 

2/o Theoretical zero-point of turbulent velocity distribu- 
tion 12 ' (16) 

z Exponent of suspended distribution 17 (27) 

a, etc & a' Angles 4 

0 A logarithmic function 36 (51) 

0x A logarithmic function 36 (49) 

y An angle 5 

8 The thickness of the laminar sublayer., 8 (5) 

A The apparent roughness diameter 8 (3) 

h' The laminar sublayer thickness for u» 10 (10) 

17 Variability factor of lift 36 (46) 

jjo Root-mean-square value of tj 37 

if* v measures in ij 0 values 37 

X Single step of bedload measured in diameters 34 

p Kinematic viscosity _ 8 

£ "Hiding factor" of grains in a mixture 36 (51) 

si Density of the fluid 10 

f a Density of the solids- 10 

r A shear stress 16 

to The shear stress at the wall or bed 8 (4) 

r y Shear stress at a distance y from the wall or bed 16 

* Intensity of transport 34 (42) 

** Intensity of transport for individual grain size 34 (41) 

* Intensity of shear on particle 36 (49) 

9* Intensity of shear on representative particle 10 (11^ 

Intensity of shear for individual grain size 37 (54) 
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Work Charts 



Figure 1— Function ^=0.216 j0r^f f A (~^)' dY ™ terms of A 



for various values of z —zrf-± 

M-\ ri /i y\z 

Figure 2 —Function 7 2 =0.216 ^_ A y J A lo g« (y) \y) ^ in terms 

of A for various values of z = l .' + 

Figure 4. — Correction x in the logarithmic friction formula in terms 

of k s /8. 

Figure 5. — Friction u" due to channel irregularities. 

Figure 6. — Settling velocity V„ for various sizes of quartz grains 

according to Rubey. 

Figure 7. — Pressure reduction in sublayer. 

Figure 8. — Pressure correction in the transition to a smooth bed. 



Figure 9. — curve compared with measured points for uniform 

sediment. 



Figure 10. — < £*- , J'* curve. 
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Figure 9. ^ curve 

compared with measured points 
(or uniform sediment. 
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